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ABSTRACT 

The ACIS instrument on board the Chandra X-ray Observatory has been used to carry out the first 
systematic study of low-mass X-ray binaries (LMXBs) in M87, the giant elliptical galaxy near the 
dynamical center of the Virgo Cluster. These images — having a total exposure time of 154 ks — are 
the deepest X-ray observations yet obtained of M87. We identify 174 X-ray point-sources, of which ~ 
150 are likely LMXBs. This LMXB catalog is combined with deep F475W and F850LP images taken 
with ACS on HST (as part of the ACS Virgo Cluster Survey) to examine the connection between 
LMXBs and globular clusters in M87. Of the 1688 globular clusters in our catalog, fx — 3.6±0.5% 
contain a LMXB. Dividing the globular cluster sample by metallicity, we find that the metal- rich 
clusters are 3±1 times more likely to harbor a LMXB than their metal-poor counterparts. In agreement 
with previous findings for other galaxies based on smaller LMXB samples, we find the efficiency of 
LMXB formation to scale with both cluster metallicity, Z, and luminosity, in the sense that brighter, 
more metal-rich clusters are more likely to contain a LMXB. For the first time, however, we are 
able to demonstrate that the probability, px, that a given cluster will contain a LMXB depends 
sensitively on the dynamical properties of the host cluster. Specifically, we use the HST images to 
measure the half-light radius, concentration index and central density, p$, for each globular, and 
define a parameter, T, which is related to the tidal capture and binary-neutron star exchange rate. 
Our preferred form for px is then px oc rp^ ' 42±011 (Z/Zo) - 33±01 . We argue that if the form of 
px is determined by dynamical processes, then the observed metallicity dependence is a consequence 
of an increased number of neutron stars per unit mass in metal-rich globular clusters. Finally, we 
present a critical examination of the LMXB luminosity function in M87 and re-examine the published 
LMXB luminosity functions for M49 and NGC 4697. We find no compelling evidence for a break 
in the luminosity distribution of resolved X-ray point sources in any of these galaxies. Instead, the 
LMXB luminosity function in all three galaxies is well described by a power law with an upper cutoff 
at L x ~ 10 39 erg s -1 . 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: individual (M87) — galaxies: star 
clusters — globular clusters: general — X-rays: binaries 



1. INTRODUCTION 

Normal elliptical galaxies have long been known to har- 
bor two major components of X-ray emission — a soft 
component due to emission from diffuse gas, and a harder 
one arising from a population of low-mass X-ray bina- 
ries (LMXBs). The existence of the latter component 
was first inferred from the spectral hardening of ellipti- 
cal galaxies with progressively smaller X-ray to optical 
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luminosities, a trend reminiscent of late-type galaxies in 
which a portion of the X-ray emission could be identi- 
fied directly with a population of accreting binary stars 
(Kim, Fabbiano & Trinchieri 1992). With the launch 
of Chandra, the harder X-Ray component has been par- 
tially resolved into point sources associated with LMXBs 
for an ever-increasing number of early-type galaxies (e.g., 
Sarazin, Irwin & Bregman 2001, 2002; Angelini, Loewen- 
stein & Mushotzky 2001; Blanton, Sarazin & Irwin 2001; 
Finoguenov & Jones 2002; Kundu, Maccarone & Zepf 
2002; Maccarone, Kundu & Zepf 2003; Kundu et al. 
2003; Jeltema et al. 2003; Kim & Fabbiano 2003; Ir- 
win, Athey & Bregman 2003; Sivakoff, Sarazin & Irwin 
2003). 

In the Milky Way, it was realized soon after the discov- 
ery of X-ray emission from a handful of Galactic globu- 
lar clusters (GCs) (Giacconi et al. 1974; Forman et al. 
1978) that the number of X-ray sources per unit mass 
is several hundred times higher in GCs than in the halo 
field (Katz 1975; Clark 1975). The population of X-ray 
point sources associated with GCs was later found to 
be a mixture of dim (Lx < 10 34 5 erg s _1 ) and bright 



(10 3b £ L x £ 10 38 erg s" 1 ) sources (Hertz & Grindlay 
1983). The fainter sources seem to be composed of many 
different kinds of objects (see, e.g., Verbunt 2001), while 
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the bright sources are believed to be accreting neutron 
stars, and are thus classified as LMXBs. Since GCs con- 
tain only <^ 0.1% of the Galaxy's stars, but <~ 10% of its 
LMXBs, it is clear that GCs are efficient sites of LMXB 
formation. According to current thinking, the overabun- 
dance of LMXBs in GCs is a direct consequence of their 
high central densities. In such environments, the rates 
of tidal capture of neutron stars and single-binary ex- 
change interactions — the two principal mechanisms by 
which LMXB progenitors are thought to be produced — 
are greatly enhanced relative to the field (Clark 1975; 
Fabian, Pringle & Rees 1975; Hills 1976). 

Chandra's ability to resolve LMXBs in nearby galax- 
ies allows us to examine the connection between LMXBs 
and GCs in new and different environments. Such stud- 
ies may yield valuable information on the processes by 
which LMXBs form and evolve. Moreover, what is lost 
in the detailed description of individual LMXBs and 
GCs beyond the Local Group is counter-balanced by 
the potentially dramatic gains in sample size (for exam- 
ple, the total number of bright X-ray sources belonging 
to the Galactic GC system is limited to just thirteen 
objects; Verbunt 2001, White & Angelini 2001). M87 
(NGC 4486), the giant elliptical galaxy near the dynami- 
cal center of the Virgo Cluster, has the richest GC system 
in the local supercluster with a total of 13,450±950 GCs 
(McLaughlin, Harris & Hanes 1994). It is also one of 
the most thoroughly studied GC systems, with a wealth 
of spectroscopic and photometric information available 
on its metallicity distribution, spatial structure, luminos- 
ity function, age distribution, and dynamical properties 
(McLaughlin, Harris & Hanes 1994; Cohen, Blakeslee & 
Ryzhov 1998; Harris et al. 1998; Kundu et al. 1999; 
Hanes et al. 2001; Cote et al. 2001; Kissler-Patig, Brodie 
& Minniti 2002; Jordan et al. 2002). Yet nothing is 
known about the LMXB population in M87 or its rela- 
tion to the M87 GC system, partly due to the difficulty 
of detecting individual X-ray point sources superimposed 
on a bright, and spatially varying diffuse background. 
In this paper, we combine deep X-ray observations from 
Chandra with optical 3475 and z 8 50 imaging from HST 
to characterize the LMXB population in M87, and to 
examine its connection to the underlying GC system. 

In what follows, we assume a distance to M87 of 
D = 16 Mpc (Tonry et al. 2001), an effective radius of 
R e = 96" (de Vaucouleurs & Nieto 1978), and a Galactic 
column density of Nh = 2.5 x 10 20 cm~ 2 (Stark et al. 
1992). Before starting, we set some notation and con- 
ventions. We will make repeated use of the Kolmogorov- 
Smirnov (KS) test and the Wilcoxon rank sum test 
(Wilcoxon 1945; Mann & Whitney 1947). The KS test 
tests the hypothesis that two samples are drawn from 
the same parent continuous distribution (two-sample KS 
test) or that a given sample was drawn from a specified 
continuous distribution (one-sample KS test), whereas 
the Wilcoxon rank sum tests the hypothesis that the lo- 
cation of the samples' parent distributions is the same. 
The results of these tests will be often reported by their 
returned p- values, which give the probability of obtaining 
the observed statistic under the null hypothesis. When 
talking about probability density functions the symbol 
~ should be understood in the sense of "distributes as" 
rather than its usual sense, and we will not write explic- 



itly in those distributions the normalization constants. 

2. OBSERVATIONS AND DATA REDUCTIONS 
2.1. X-ray catalog 

M87 (=NGC 4486) was observed with the Chandra 
Advanced CCD Imaging Spectrometer (ACIS) for 121 
ks on 5-6 July 2002. In what follows, we use only the S3 
chip data. The data were processed following the CIAO 
data reduction threads, including a correction for charge 
transfer inefficiency (CTI; Townsley et al. 2000). Addi- 
tionally, we used 38 ks of archival ACIS observations of 
M87 taken on 29 July 2000 (PI: A.S. Wilson). These data 
were processed in similar fashion to the July 2002 data, 
except that no CTI correction was possible because the 
data were telemetered in Graded mode. All reductions 
were carried out with CIAO 2.3 coupled with CALDB 
2.21. In order to combine the events files into a single 
image for point source detection, we obtained relative 
offsets by matching the celestial coordinates of two X- 
ray point sources 9 . The relative offset was w 0'.'5. After 
excluding the provisional list of point sources identified 
with SExtractor (Bertin & Arnouts 1996), the M87 jet, 
and the central regions of the galaxy, we created a light 
curve binned in 50 s time intervals in order to clean back- 
ground flares by rejecting points deviating by more than 
3(J from the quiescent mean. The total exposure time of 
the co-added image, excluding four flares totaling «2.5 
ksec, was 154 ks. 

Detection was performed using a two stage process. 
First, a wavelet filtering of the image was performed 
to keep only those objects with a characteristic struc- 
ture <; 8 pixels and a significance greater than 2 x 10~ 5 
(~ 4er) using the MR/1 package (Starck, Murtagh & Bi- 
jaoui 1998). Object detection was then performed on 
this filtered image using SExtractor. Note that this de- 
tection procedure was also used by Valtchanov, Pierre & 
Gastaud (2001) who found it to be the most effective for 
XMM-Newton data. The inner regions of M87 exhibit a 
wealth of structure in its diffuse emission, with numerous 
bubbles and arcs that cause spurious point source detec- 
tions. These regions, along with the jet and the central 
source, were masked. Source regions were defined to be 
ellipses with semi-major and semi-minor axes equal to 
twice the parameter values returned by SExtractor. We 
also used wavdetect (Freeman et al. 2002) to carry out 
the source detection, producing a similar point source 
catalog in the process, but a visual inspection of the in- 
dividual detections suggested that the adopted method 
produced superior determinations of the source centroids 
(which are of critical importance when matching to the 
optical sources). Indeed, using positions obtained with 
wavdetect the rms of the difference in celestial coordi- 
nates between the X-ray and optical sources increased 
by - 35%. 

The total number of detected point sources in the S3 
chip was 174, of which <~ 20 are expected to be back- 
ground contaminants such as AGN (Mushotzky et al. 
2000; Giacconi et al. 2001). A background region was 
defined for each source by taking an annulus centered on 
the source. Since the background varies on small spatial 
scales in the inner regions of the galaxy, the inner radii 

9 CXOU J123047.1+122415 and CXOU 123044.6+122140 
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of the annuli were taken to be 

n = a ma j for r < 0.44i? e , . 

r t = l.5a maj for r > O.UR e [ > 

where a ma j is the semi-major axis of the source extrac- 
tion region. For all objects, the outer radius of the back- 
ground annulus was calculated so that the total back- 
ground area was five times that of the source extraction 
region. For a few objects, the background annulus de- 
fined in this way included a second point source. In such 
cases, the annulus was masked over the angular region 
containing the adjacent source, and the outer radius ad- 
justed so that the area of the background region was five 
times that of the source extraction region. 

2.2. Optical catalog 

M87 was observed as part of the ACS Virgo Cluster 
Survey (Cote et al. 2004, in preparation) on 19 January 
2003. The full dataset consists of two 375 s exposures 
in the F475W (3475) band, two 560 s exposures in the 
F850LP (z85o) band, and a single 90 s F850LP exposure. 
A detailed account of the reduction procedures for the 
survey will be presented elsewhere (Jordan et al. 2004, 
in preparation) so we give only a brief summary here. 

After determining the small shifts between the images, 
they were combined using the PYRAF task multidrizzle 
(Koekemoer et al. 2002). Detection images are built by 
subtracting a model of the galaxy and object detection 
on these images was performed independently for each 
band using SExtractor (Bertin & Arnouts 1996) with a 
threshold of 5 connected pixels at a level of 1.5 a. The 
celestial coordinates of the detections were matched with 
a O'.'l matching radius; sources detected in just a single 
filter were discarded. 

GCs at the distance of M87 are slightly resolved with 
ACS. This opens the possibility of modeling directly the 
two-dimensional light distribution of the GCs. A code 
has been developed (Jordan & Cote 2004, in preparation) 
to measure total magnitude, half-light radius, Vh, and 
concentration index, c, for each GC by fitting the two- 
dimensional ACS surface brightness profiles with the con- 
volution of the instrumental point spread function (PSF) 
with isotropic, single-mass King (1966) models. These 
models are well known to provide an excellent represen- 
tation of the surface brightness profiles of most Galactic 
GCs. DAOPHOT II (Stetson 1987; 1993) was used to 
derive PSFs that varied quadratically with CCD posi- 
tion, using archival observations of moderately crowded 
fields in the outskirts of the Galactic GC 47 Tucanae. 
Instrumental magnitudes were converted to the AB sys- 
tem using zeropoints of 26.07 mag for F475W and 24.86 
mag for F850LP (Sirianni et al. 2004). A correction for 
foreground extinction was performed using the redden- 
ing curves of Cardelli, Clayton & Mathis (1989), with a 
value of E(B-V) = 0.023 taken from the DIRBE maps of 
Schlegel, Finkbeiner & Davis (1998). Hereafter, we de- 
note the F475W and F850LP filters by the corresponding 
filters (3475 and zsso, respectively) in the Sloan Digital 
Sky Survey. 

In order to define a set of bonafide GCs using this opti- 
cal catalog, several additional selection criteria were im- 
posed: (1) a color in the range 0.4 < (3475 — ^sso) < 1-9, 
as measured from both SExtractor and the King model 
fitting program; (2) a half-light radius in the range 



0.5 < r/i(pc) < 15.5, an interval that encompasses almost 
all GCs in the Milky Way (Harris 1996); and (3) agree- 
ment between the measurements in the 3475 and z^5o 
bandpasses at the 4<r level. Additionally, we removed 
from the optical catalog 106 GC candidates which fall 
in the regions that were masked in the X-ray image (see 
§ I2.1|) . A total of 1688 GC candidates, spanning nearly 
six magnitudes in brightness, met these selection crite- 
ria. When matching to the X-ray point-source catalog, 
the full list of optical point sources has been used, since 
it is likely that some fraction of the X-ray point sources 
will be unrelated to GCs. 

2.3. Matching 

The high density of GCs in M87 requires care to be 
taken in matching the X-ray and optical catalogs. As an 
additional complication, the factor of ten difference in 
pixel size between ACIS and ACS makes any uncertainty 
in the X-ray coordinates translate into many ACS pixels. 

The quality of the absolute pointing was first ver- 
ified by comparing the coordinates of the galaxy nu- 
cleus, in both the X-ray and optical images, with the 
VLBI position of the M87 nucleus obtained by Ma et al. 
(1998): a(J2000)= 12 /l 30 m 49. s 423381 and 5(J2000)= 
12°23'28."0434. For both the optical and X-ray images, 
the coordinates of the brightest nucleus pixel differ by 
less than ~ 0'.'3 from the VLBI coordinates. We are thus 
confident that both datasets have good absolute point- 
ing. An initial matching without any adjustment be- 
tween the X-ray and optical sources further confirmed 
the compatibility of the coordinates, and an analysis of 
the residuals showed no statistically significant trends as 
a function of position in the chip. 

Given that any offset is small and there is no appre- 
ciable rotation between the two coordinate systems, we 
adopted the following iterative scheme to obtain constant 
offsets in each coordinate, Aa and AS. First, all sources 
within a radius of 0'.'3 were first matched, and a biweight 
(Beers, Flynn & Gebhardt 1990) estimate of the offsets 
calculated. In the next iteration, the matching radius 
was then incremented to 0'.'5, improved biweight esti- 
mates of the offsets calculated, and the list of matched 
objects updated accordingly. This process was repeated 
until the list of matched objects stabilized. 

The adopted offsets (in the sense optical minus X-ray) 
are Aa = -O'.'l and AS = +0'.'15. The rms deviation 
of the residuals are (X'14 and 0'.'13, respectively. Fig- 
ure H shows the central 5' x 5' of the Chandra image 
with the ACS field of view overlaid. X-ray point sources 
are indicated by green ellipses; those sources that coin- 
cide with optical GC candidates are marked with white 
squares. The final list contains 62 optical sources of any 
sort that are matched to an X-ray source; 60 of these 
optical sources are probable GCs based on the criteria 
described in W2.2\ Data for all X-ray point sources are 
given in Tabled which records the source identification 
number, coordinates, count rate, luminosity and hard- 
ness ratios (see below) . Column 7 gives a flag to indicate 
whether the X-ray source falls within the ACS field of 
view, while comments on the classification of the various 
optical sources are given in the final column. 

Two sources in particular merit attention (see Fig- 
ure |2J). The extraction ellipse of one X-ray source, 
CXOU J123047. 1+122415 in Table □ encloses three op- 
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Fig. 2.— (Left Panel) Portion of the ACS F850LP im- 
age in which the extraction region for the X-ray source CXOU 
J123047. 1+122415 is marked by an ellipse. The X-ray source con- 
tains three GC candidates, none of which is within a matching 
radius of the X-ray emission centroid. (Right Panel) Portion of the 
ACS F850LP image in which the extraction region for the X-ray 
source CXOU 123046.7+122402 is marked by an ellipse. The X-ray 
source contains two GC candidates within one matching radius of 
its centroid. 



Fig. 1 — Co-added Chandra/ACIS image of M87 with the ACS 
field of view overlaid (rhomboids). The X-ray point source detec- 
tions are indicated by the green ellipses. White squares indicate 
the 60 X-ray sources that coincide with globular cluster candidates. 
North is up and East is to the left in this image, which measures 
5' x 5'. 



tical sources. In X-rays, the source is extended in a way 
that is consistent with being a blend of multiple sources. 
The centroid of the X-ray emission lies approximately 
halfway between the brightest optical detections, yet the 
three optical sources lie outside the nominal matching 
radius and so do not make it into the final catalog. We 
consider this to be a match when calculating the frac- 
tion of X-ray-optical matches by assuming that two of 
the GC matches hold an LMXB, but discard this source 
from the subsequent analysis. A second X-ray source 
(object CXOU J123046. 7+122402 in Table[TJ has two op- 
tical candidates within the matching radius. Given this 
ambiguity, we consider this to be a match for the pur- 
poses of estimating the overall frequency of LMXB-GC 
associations, but do not include this source in any other 
aspect of the analysis. All the candidate GCs in these 
two sources are metal-rich, so no ambiguity is introduced 
when calculating the frequencies for the metal-rich and 
metal-poor groups. Removing the two optical matches to 
CXOU J123046. 7+122402 leaves 58 X-ray point source 
matches that will be used in the analysis. 

Given the high GC surface density in M87, it is of in- 
terest to know the number of chance matches that might 
occur between the X-ray and optical datasets. We have 
estimated the number of such false matches by rotating 
the X-ray source coordinates about galaxy's nucleus and 
calculating the total number of matches at each rota- 
tion angle. This exercise produced an average of four 
matches at each angle, so we conclude that the number 
of chance associations in our sample is small. According 
to the background source counts of Giacconi et al. (2001) 
and Mushotzky et al. (2000), we expect ~ 2 background 
sources within the ACS field of view. This is comparable 
to the number (2) of X-ray sources in our sample that 
match an optical source that is not a probable GC, so 
we believe that our sample has very little contamination 
from background sources and false matches. 



2.4. Variability 

We performed a simple search for variability among the 
23 X-ray sources with Lx > 10 38 erg s^ 1 that coincide 
with a GC candidate. Using the positions from the com- 
bined observations, we measured the fluxes for the 2000 
and 2002 datasets. The distribution of the flux differ- 
ences, when divided by the expected uncertainty, reveals 
that seven sources show significant variability when com- 
pared with a normal distribution. Thus, the incidence of 
X-ray transients in the M87 GC/LMXB population ap- 
pears roughly consistent with that of the Galaxy, where 
roughly half of all sources are known to exhibit time vari- 
ability (e.g., Verbunt et al. 1995). 

3. SPECTRAL ANALYSIS 

All X-ray sources with galactocentric radii in the range 
0.44i? e < r < 2R e , were extracted and summed into 
a composite spectrum with the CIAO routine acisspec, 
which computes a weighted redistribution matrix file 
and ancillary response file (ARF) appropriate for point 
sources distributed over the ACIS detector. The ARF 
files were corrected for the degradation of the A CIS 
quantum efficiency using the CIAO tool apply -acisabs, 
which applies the ACISABS absorption profile (Chartas 
& Getman 2002) 10 to the original ARF file. The extrac- 
tion was performed independently for the 2000 and 2002 
datasets, and the source spectra were regrouped so that 
each energy channel contained a minimum of 50 pho- 
tons, prior to background subtraction. In Figure wc 
show summed spectra from the 2002 dataset, for both 
the source regions (i.e., object plus background) and for 
the background regions alone. 

The composite, background-subtracted source spectra 
for the 2000 and 2002 datasets were fit simultaneously 
with XSPEC vll.2.0 (Arnaud 1996). The spectral en- 
ergy distribution is assumed to have a single photon 
power-law behaviour oc E~ K . The Galactic hydrogen 
column density is held fixed at Njj = 2.5 x 10 20 cm -2 
in the fit; channels with energy less then 0.7 keV and 
greater than 4 keV were excluded from the fit. The spec- 
tra, best-fit model and residuals are shown in Figure 21 
The best-fit power-law exponent is n = 1.641q'q46 (90% 

10 See http: //www. astro .psu. edu/users/chartas/xcontdir/xcont .html 
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Fig. 3. — Summed spectra for the background regions (red 
crosses) and for source regions prior to background subtraction 
(black crosses). The spectra are based entirely on data acquired in 
2002. 





Fig. 5. — Left: Hardness ratios, H21 and H31, for the full sample 
of 174 X-ray point sources. The size of the symbol for each source 
is proportional to its luminosity in the range 0.3-10 keV. The mean 
hardness ratios are ((H21), (H31)) = (0.07, -0.21), as indicated by 
the cross. Right: Hardness ratios, H21 and H31, for all sources 
with Lx > 1-5 X 10 38 erg s — 1 . Error bars give la uncertainties in 
the ratios. The line shows predicted ratios for a power law with a 
hydrogen column density equal to Njj = 2.5 X 10 20 cm -2 (Stark 
et al. 1992); from top to bottom, the squares correspond to the 
predictions for power law exponents of 0, 1, 2 and 3. 
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Fig. 4. — Summed spectra for all X-ray sources, with the 2000 
data shown in black and the 2002 data in red. The curves indicate 
the best-fit power-law model, with index re = 1.64, obtained by 
fitting simultaneously to both data sets. The lower panel shows 
the residuals from this best-fit model. 

confidence uncertainties), with a reduced chi-squared of 
xt = 1.032 with 218 degrees of freedom. Irwin et al. 
(2003) analyzed the composite spectra of point sources 
within 3i? e for a sample of 15 nearby galaxies, finding 
power-law exponents in the range 1.45 < k < 1.9. Thus, 
our measured power-law exponent for the LMXBs in M87 
is typical of those found in other nearby galaxies. This 
best-fit model is used to convert the observed counts to 
unabsorbed luminosities, Lx , over the range 0.3—10 keV, 
assuming that all the sources are at the distance of M87. 
The resulting conversion factor is 1.4 x 10 41 erg count -1 . 

We examined the crude spectral properties of the re- 
solved sources by calculating hardness ratios, which have 
the advantage of being measurable for even the faintest 
sources. Counts were calculated for three distinct energy 
bands: a soft (0.3 — 1 keV) band denoted by S, a medium 
band (1-2 keV) denoted by M, and a hard (2 - 10 keV) 
band denoted by H. Following Sarazin et al. (2000), the 
hardness ratios, H31 and H21, are taken to be 

H21 = (M — S)/ (M + S) r v 

H31 = (H — S)/(H + S) [Z > 

The distribution of hardness ratios is shown in Figure [SJ 



The sources occupy a diagonal swath in this plot, as is 
typical for LMXBs (Sarazin et al. 2000; Blanton et al. 
2002; Irwin et al. 2003; Jeltema et al. 2003; Sivakoff et al. 
2003). The mean location, at (0.07,-0.21), is indicated 
by the cross. It is apparent from this figure that the most 
luminous sources appear to have the softest spectra, and 
a Wilcoxon test confirms this impression, giving proba- 
bilities of 6% and 0.6%, respectively, that the H21 and 
H31 values for sources with Lx > 5 x 10 38 erg s _1 share 
the same location. Although there are only six objects 
with Lx > 5 x 10 38 erg s _1 , this result is consistent with 
that of Irwin et al. (2003), who noted that sources with 
Lx > 10 39 erg s~ x appear to be significantly softer. As 
they remark, this trend might be a reflection of inverse 
dependence between the emitted flux and spectral state 
exhibited by candidate black hole X-ray binaries in the 
Milky Way {e.g. Tanaka & Lewin 1995; Nowak 1995). 
Finally, we note that there are two sources in the upper 
right corner of Figure [S] in which both H21 and H31 are 
equal to 1.0; given the hardness of these spectra, we sus- 
pect that these sources may be strongly absorbed AGN. 

4. GLOBAL PROPERTIES OF THE LMXB POPULATION 

We now turn our attention to the observed properties 
of the LMXB population as a whole (i.e., spatial struc- 
ture, luminosity function, and suitability as a distance 
indicator). Our ultimate aim is to understand the nature 
of the connection between LMXBs and GCs in M87, and 
to examine the broader implications for LMXB forma- 
tion. 

4.1. Radial and Azimuthal Structure 

How does the radial distribution of LMXBs in M87 
compare to that of its GCs and the underlying galaxy 
light? In comparing the various profiles, we restrict our- 
selves to X-ray point sources having Lx > 7 x 10 37 ergs 
s _1 , and to GC candidates that do not fall within the 
regions masked during the X-ray point source detection 
procedure and having z 850 < 22.8 mag (a total of 867 
objects) in order to guard against incompleteness effects. 
The heavy solid curve in Figure El shows the resulting cu- 
mulative distribution for the M87 GC system, calculated 
directly from the GC catalog described in H2.2I 



6 



LOW MASS X-RAY BINARIES AND GLOBULAR CLUSTERS IN M87 



In principle, it should also be possible to measure the 
profile of the galaxy itself from our ACS images. How- 
ever, such an approach is undermined by the limited 
areal coverage of the ACS field and the fact that our 
ACS images, which are centered on the galaxy's nucleus, 
provide limited constraints on the background surface 
brightness. Given these problems, we estimate the cu- 
mulative light distribution within the ACS field, S(r), 
by using the wide-field surface photometry of Caon, Ca- 
paccioli & Rampazzo (1990). For an annulus centered 
on the galaxy, the fractional area falling within the ACS 
field (excluding those regions masked during X-ray point 
source detection; see below) is f a (r). The cumulative 
distribution of galaxy light is then 

S(r) = 27rr/ a (r)10-°- 4 ^M (3) 

where ^i(r) refers to the -B-band surface brightness profile 
of Caon et al. (1990). Note that the implicit assumption 
of circular symmetry in Equation[3]is quite reasonable for 
M87, which has a luminosity-weighted mean cllipticity of 
(e) ~ 0.05 inside r ~ 2.'75 — the maximum galactocen- 
tric radius of our ACS field. 

The thin solid curve in Figure shows the cumula- 
tive profile of the galaxy light within the ACS field. A 
KS test confirms the well known result that the GC sys- 
tem of M87 has a shallower profile than the galaxy itself 
(e.g., Grillmair et al. 1986; Harris 1986). Also shown in 
Figure El are the cumulative distributions for two LMXB 
subsamples: the dotted curve shows the distribution for 
those sources that are associated with GC candidates, 
while the dashed curve shows the distribution for the re- 
maining X-ray sources. In both cases, we plot only those 
X-ray sources that fall within the ACS field. Since both 
samples are subject to the same selection effects, it is 
straightforward to compare these distributions directly. 
A two sample KS test accepts the hypothesis that they 
were drawn from the same parent sample. 

Comparing the X-ray point source samples to the 
galaxy and GC profiles is more difficult. Our Chandra 
image reveals the inner ~ 40" of M87 to have a remark- 
ably complex structure in diffuse emission (Figure ^ see 
also Figure 1 of Young, Wilson & Mundell 2002 and 
Sparks et al. 2004). Because of this complexity, it was 
necessary to mask several problematic regions prior to 
object detection (see limiting the region over which 
the various profiles can be compared. Perhaps as a con- 
sequence, a one-sample KS test accepts the hypothesis 
that both X-ray samples (i.e., those point sources with, 
and without, an associated GC candidate) were drawn 
from the same parent distribution as the galaxy light; 
moreover, a two-sample KS test accepts the hypothesis 
that they were drawn from the same parent distribution 
as the GC candidates. Stronger conclusions will require 
an expanded census of LMXBs but, given the brightness 
and complexity of the diffuse X-ray emission in the in- 
ner regions of M87, the requisite observations will prove 
extremely challenging. 

We also explored the azimuthal distribution of the X- 
ray point source samples. In Figure[3we show the cumu- 
lative distribution function for the X-ray point sources 
that are associated with a GC and those that are not. 
We also show the azimuthal distribution of the full GC 
sample (excluding regions that were masked in the X- 
ray point source detection). It is clear from the figure 
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Fig. 6. — Comparison for the radial distribution of galaxy light, 
globular clusters and LMXBs in M87. The heavy solid curve shows 
the cumulative distribution of 867 globular cluster candidates in 
the ACS field that are brighter than zg50 = 22.8. The thin solid 
curve shows distribution of galaxy light within our ACS field. The 
dotted curve shows the cumulative distribution of the 32 LMXBs 
in our ACS field that are brighter than Lx = 7 X 10 37 erg s _1 and 
coincide with a globular cluster candidate. The dashed curve show 
the cumulative distribution of the remaining sample of 15 LMXBs 
in the ACS field that are brighter than Lx = 7 X 10 37 erg s _1 but 
have no GC counterpart. 



that the X-ray point sources associated with a GC fol- 
low the distribution of the full GC sample, and this is 
confirmed by a KS test. The X-ray point sources not 
associated with a GC seem somewhat deviant, but a KS 
test accepts the hypothesis that the sample was drawn 
from the same distribution as the full GC sample and the 
X-ray point sources associated with a GC with p-values 
°f Pks — 0.19 and pks = 0.1 respectively. This exer- 
cise reveals that the spatial distribution of X-ray point 
sources associated with a GC is representative of the full 
GC sample. 

4.2. Luminosity Function 

The luminosity function of LMXBs is of considerable 
interest, both as a rare constraint on the mass distri- 
bution of accreting sources in external galaxies, and as 
a potential distance indicator. Working with a sample 
of ps 80 LMXBs in NGC 4697, Sarazin et al. (2001) 
showed that their cumulative luminosity function was 
well described by a broken power-law, with a "break" 
at Lb ps 3.2 x 10 38 ergs s . Since this is close to the Ed- 
dington luminosity for spherical hydrogen accretion onto 
the surface of a 1.4M Q neutron star (e.g. Shapiro & 
Teukolsky 1983), Sarazin et al. (2001) drew attention to 
the possibility of using this feature as a standard candle 
in distance estimation. Indeed, the use of a character- 
istic luminosity in accreting neutron stars as a distance 
indicator dates to early studies of Galactic X-ray sources 
(e.g., Margon & Ostriker 1973; van Paradijs 1978). 

4.2.1. Representation as a Broken Power-Law 

In Figure |H| we show the cumulative luminosity func- 
tion, log N(Lx), of X-ray sources in M87. Here, N(Lx) 
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Fig. 7. — Comparison for the azimuthal distribution of globular 
clusters and LMXBs in M87. The solid curve shows the cumulative 
distribution of all bona fide GC candidates. The dashed curve 
shows the cumulative distribution of the X-ray point sources that 
coincide with a GC and the dotted curve shows the cumulative 
distribution of the X-ray point dources that do not coincide with 
a GC. 
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Fig. 8. — Observed log N(Lx), where N(L X ) is the number of X- 
ray sources with L > Lx ■ The dashed line shows the distribution 
for all sources within the S3 chip; the sources within the ACS 
field and associated with globular cluster are represented by the 
solid line; those sources within the ACS field but not associated 
with a globular clusters are shown by the dot-dashed line. The 
smooth lines are the best-fit power-law luminosity distributions 
of the form / oc L~^ J with an upper cutoff such that / = for 
Lx > 10 39 erg s — 1 . For the full sample of all sources, and for the 
subsample of sources that not associated with a globular cluster, 
the background contamination has been removed as in Giacconi 
et al. (2001). 



is the number of objects with luminosities in excess of 
Lx- The dashed line shows this distribution for the com- 
plete sample of sources — a total of 174 objects. The 
two lower distributions show the distributions for those 
sources that fall within the ACS field: the solid line 
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Fig. 9. — Same as Figure 151 except for sources in M49 and 
NGC 4697. The dashed histogram shows the distribution of X-ray 
sources in M49 from the S3 chip of Kundu et al. (2002). The 
corresponding distribution for sources in NGC 4697 is given by the 
solid histogram (Sarazin et al. 2001). The dot-dashed histogram 
shows the distribution for sources in M49 that are associated with a 
globular cluster (Kundu et al. 2002). The smooth curves are best- 
fit power-law luminosity distributions of the form / oc L~^' with 
an upper cutoff such that / = for Lx > 10 39 erg s — 1 . For the 
complete samples in both galaxies, the background contamination 
has been removed as in Giacconi et al. (2001). 



indicates those sources that coincide with a GC, while 
the dot-dashed line refers to sources with no associated 
GC. Aside from normalization, the three distributions 
appear remarkably similar; the abrupt flattening below 
Lx ~ 4 x 10 37 erg s _1 is probably due to incomplete- 
ness. When fitting models to the luminosity functions, 
we have corrected for background contamination using 
the background counts of Giacconi et al. (2001) for the 
complete sample and the subsample of objects that are 
not associated with GCs. 

To facilitate comparison with previous work, we have 
fitted these three distributions with broken power-laws 
of the form 



N(L X ) <x 
N(L X ) <x L° 2 



r oi 



for 
for 



Lx < Lb 
L x >L b 



(4) 



To guard against incompleteness effects, we consider 
only those source brighter than 5 x 10 37 erg s . The 
resulting values for a±, a-i and Lb are given in Ta- 
ble [21 Note that the break luminosities found here, 
Lb ~ (2 — 3) x 10 38 erg s _1 , are similar to those found 
in other early-type galaxies (e.g., Sarazin et al. 2001; 
Finoguenov & Jones 2001; Kundu et al. 2002). Although 
the observed luminosity distribution is well described by 
this particular choice of parameterization, the data do 
not require a broken power-law. As we now show, the 
data are equally well represented by a single power-law 
distribution with an upper cutoff in luminosity. 

4.2.2. Representation as a Truncated Power-Law 

For each of the three samples of observed X-Ray lumi- 
nosities, whose corresponding cumulative distributions 
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Fig. 10. — The smooth solid line shows the expected behaviour 
of log N (> L) for a sample of size 100 drawn from a distribution of 
luminosities given by L -2 with an upper cut at L = 10 39 erg s — 1 . 
The rest of the lines show the observed behaviour of 5 random 
realizations of size 100 from that distribution. 



are shown in Figure we have fit single power-law dis- 
tributions of the form L x for L m in < Lx < L max , taking 
L m ax = 10 39 erg s _1 as the upper cutoff. As before, we 
adopt L m in — 5 x 10 37 erg s _1 to guard against incom- 
pleteness at the faint end. The corresponding best-fit cu- 
mulative distributions are shown as the smooth curves. 
In each case, a one-sample KS test reveals that the sam- 
ples are consistent with their being drawn from a single 
power law with 7 ~ —2.1. The best fit values for 7 are 
given in Table ET 1 - 

Is this true of the LMXB populations in other galax- 
ies? In Figure |5| we show the cumulative luminosity 
functions of all LMXBs in NGC 4697 and M49 (using 
data from Sarazin et al. 2001 and Kundu et al. 2002, 
respectively). Maximum likelihood fits to both datasets 
reveals that they are consistent with having been drawn 
from a single power-law distribution with an upper cut 
at Lx = 10 39 erg s _1 . The corresponding best-fit cu- 
mulative distributions, which have 7 = —1.63 ± 0.14 
(NGC 4697), 7 = -1.58 ± 0.12 (M49, all X-ray point 
sources) and 7 = —1.64 ± 0.25 (M49, X-ray sources as- 
sociated with a GC), are given by smooth curves in Fig- 
ure Ej] As this exercise demonstrates, it is dangerous 
to draw conclusions about the underlying distribution, 
/, from the quantity N(Lx), particularly at the high- 
luminosity end, if / is truncated above some value. Let 
us denote the cumulative distribution of / by F, so that 
log N{Lx) is equivalent modulo a constant to log(l — F). 
The slope s of this function is given by s = — //(l — F). If 
/ = for values of Lx greater than L max then, as we ap- 
proach L max , we generically expect lhri£^£ max s = —00, 

11 Note that the broken power law is fitted using N(Lx), 
whereas the truncated power law is fitted using the observed sam- 
ples of Lx, whose parent distribution we denote by /. N(Lx) 
is (1 — F) modulo a normalization constant, where F is the cu- 
mulative distribution of /. Thus, as the bulk of the data lies at 
luminosities lower than the inferred break luminosity L^, we expect 
that ai fts 7 + 1. 



producing a dip in the expected form of log(l — F). To 
give an example germane to the present discussion, let us 
assume that / = L x 2 with / = for Lx > 10 39 erg s _1 . 
The corresponding distribution N(Lx) is shown as the 
smooth curve in Figure 1101 For comparison, the his- 
tograms in this figure show five simulated datasets, each 
consisting of one hundred objects. Because the parent 
distribution is cut beyond a maximum Lx some simu- 
lated datasets have an apparent break, even though the 
parent distribution has no characteristic scale to distin- 
guish the two regimes. 

Based on the evidence presented above, we conclude 
that there is no compelling evidence for two funda- 
mentally different accretor populations in M87, M49 or 
NGC 4697, and that a single power-law distribution 
(truncated above ~ 10 39 erg s _1 ) provides an adequate 
description of the LMXB populations in all three galax- 
ies. Moreover, we have shown that the apparent "breaks" 
at Lx ~ (1 — 4) x 10 38 erg s _1 may be an artifact 
of this distribution. Our conclusions are in agreement 
with those of Sivakoff et al. (2003), who found that 
the luminosity distribution of LMXBs in NGC 4365 and 
NGC 4382 could be better modeled by a power law hav- 
ing an upper cutoff at Lx ~ 10 39 erg s _1 . This also 
appears to be the case in the Milky Way: Grimm, Gil- 
fanov & Sunyaev (2001) find a truncated power law to 
be an excellent representation of the Galactic LMXB lu- 
minosity distribution. 

In retrospect, the lack of a characteristic luminosity 
scale should perhaps come as no surprise since the Ed- 
dington luminosity, Le, is usually computed under very 
particular assumptions: namely, spherical accretion of 
pure ionized hydrogen and Thomson scattering. This 
is clearly an idealized situation, and there are various 
ways in which an accreting neutron star can exceed this 
naive estimate: e.g., unusual chemical abundance, for- 
mation of a super-critical disk, radiation in the form 
of relativistic jets and the presence of strong magnetic 
fields (see, e.g., the discussion in Grimm et al. 2001 
and references therein). Still another effect that would 
blur a characteristic scale is that the observed values 
of Lx will be affected by disk obscuration and scatter- 
ing. Irwin et al. (2003) find essentially no sources with 
L x > 2 x 10 39 erg s" 1 in their study of the LMXB pop- 
ulations of 15 early- type galaxies, and our results are in 
line with theirs. Thus, it is likely that the luminosity 
function of LMXBs in most early-type galaxes is trun- 
cated at L max ~ (1 — 2) x 10 39 erg s^ 1 . This upper cutoff 
is probably not universal. Indeed, a significant number 
of very luminous sources (Lx <> 10 40 erg s _1 ) have been 
observed in some early-type galaxies, such as NGC 720 
(Jeltema et al. 2003). In NGC 720, the spatial distribu- 
tion suggests that the more luminous sources could arise 
from a younger stellar population whose formation was 
perhaps triggered by a recent merger. 

4.2.3. Environmental Dependence 

Do the LMXBs in GCs share the same luminosity func- 
tion as those which are not associated with GCs? There 
have been conflicting claims about such an environmental 
dependence in the LMXB luminosity function. In their 
study of NGC 1399, Angelini et al. (2001) found the 
LMXBs in GCs to be, on average, more luminous than 
those which are not associated with GCs. On the other 
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Fig. 11. — Normalized, cumulative luminosity functions for 
LMXBs within the ACS field. The solid curve shows the distri- 
bution of LMXBs that are associated with a GC; LMXBs that are 
not associated with GCs are indicated by the dashed curve. 



hand, Kundu et al. (2002) found the two populations in 
M49 to have indistinguishable luminosity functions, and 
Sarazin et al. (2003) reached a similar conclusion based 
on their analysis of LMXBs in four early-type galaxies. 

In Figure 1111 we show the normalized, cumulative lu- 
minosity function for LMXBs within the ACS field. The 
solid curve shows the distribution of LMXBs that co- 
incide with a GC, while the dashed curve shows those 
LMXBs that do not. The LMXBs associated with GCs 
are slightly brighter on average, with a mean luminosity 
of (L x ) = (1.4 ± 0.2) x 10 38 erg s -1 . The mean lumi- 
nosity of LMXBs that are not associated with GCs is 
(L x ) = (1.1 ± 0.2) x 10 38 erg s -1 . The difference, how- 
ever, is not significant at greater than 99% confidence, 
as KS and Wilcoxon sum rank test accept the hypothe- 
sis that the data were drawn from the same distribution 
with p- values of pks = 0.1 and Pwilcox — 0.1, respec- 
tively. Thus, at least in the case of M87, we find no 
support for the claim that the LMXBs that are associ- 
ated with GCs are significantly brighter than those which 
are not. These findings are consistent with the sugges- 
tion that essentially all LMXBs in early-type galaxies 
may have first formed in GCs (White, Sarazin & Kulka- 
rni 2002) and were subsequently ejected or dispersed into 
the general field. 

5. THE RELATION BETWEEN LOW MASS X-RAY 
BINARIES AND GLOBULAR CLUSTERS 

5.1. The Efficiency of LMXB Formation 

The most basic characterization of the probability that 
a GC harbors an LMXB is the overall probability, px, 
that a GC will host a LMXB candidate. Restricting our- 
selves to those sources in the ACS field of view, we find 
px = 0.036 ± 0.005 12 . This overall probability has been 

12 Because of incompleteness, this number is a slight overesti- 
mate, as we ar e mi ssing faint GCs that will not contribute many 
LMXBs (see § 15.31 . We can estimate the magnitude of the bias 
by assuming that the GC luminosity function <t>(m) is represented 



found to be in the range 2 — 4% in a wide variety of 
galaxy types (see, e.g., Sarazin et al. 2003 and refer- 
ences therein) , providing a very uniform characteristic of 
the connection between GCs and LMXBs in GCs. An- 
other basic quantity is the fraction of X-ray point sources 
fx-GC that reside in a GC; for M87 this is /x-gc = 0.62. 
Observed values of /x-gc vary substantially from galaxy 
to galaxy, from ~ 0.2 for sources in the central region 
of M31 (Primini, Forman & Jones 1993) to ~ 0.7 for 
NGC 1399 (Angelini et al. 2001), and it has been pro- 
posed that the available observations indicate a system- 
atic increases along the Hubble sequence from late to 
early types (Sarazin et al. 2003). We now turn to con- 
sider the dependence of px on various factors. 

5.2. Dependence on Metallicity 

In the Galaxy and M31, there is a tendency for 
LMXBs in GCs to appear preferentially in metal-rich 
GCs (Grindlay 1993; Bellazzini et al. 1995), but the lim- 
ited samples sizes (i.e., just 13 and 19 LMXBs, respec- 
tively, with L x > 10 36 erg s" 1 ; Verbunt 2001, White k, 
Angelini 2001) have precluded definite conclusions. How- 
ever, new observations have confirmed the trend with 
larger samples of LMXBs in NGC 1399 (Angelini et al. 
2001), NGC 4472 (Kundu et al. 2002) and NGC 5128 
(Minniti et al. 2004). Kundu et al. (2002) find that 
metal-rich GCs are, on average, ~ 3 times as likely to 
harbor LMXBs than their metal-poor counterparts. 

In Figure 1121 we show the (17475 — zsso) color distri- 
bution for the full sample of 1688 GCs in our ACS field, 
along with the corresponding distribution for the 58 GCs 
which coincide with an X-ray point source. It is appar- 
ent that LMXBs show a marked preference for metal-rich 
GCs, and this is confirmed with two different statistical 
tests: a two sample KS test rejects the hypothesis that 
the distributions were drawn from the same distribution 
with a p-value of pks = 3 x 10~ 4 , while a Wilcoxon 
rank sum test rejects the hypothesis that the parent 
distributions have the same location with a p-value of 

Pwilcox — 2 X 10 

To divide the color distribution into metal-rich and 
metal-poor subsamples, we use the KMM algorithm 
(Ashman, Bird & Zepf 1994). This algorithm performs 
the division based on an a posteriori likelihood using 
two Gaussians to model the color distribution. KMM 
gives a dividing color of (5475 — ^85o) = 1-121, which 
we adopt as the division between the metal-rich and 
metal-poor subpopulations in M87. With the subpop- 
ulations defined in this way, we find a probability of 
Px,mr = 5.1 ± 0.7% that a given metal-rich cluster will 
contain a LMXB; this should be compared to the value 
of Px,mp — I-? ^ 0.5% found for the metal-poor GCs. 13 
We conclude that metal-rich GCs are 3 ± 1 times more 
likely to contain LMXBs than metal-poor GCs. 

by a Gaussian with a = 1.4 and turnover magnitude My = —7.4 
(Harris 2001). Using this form, the luminosity function <fi of GCs 
associated with an LMXB will satisfy <f>(m) ~ L°- S9 $(m) f8 ET51 . 
where L is the GC luminosity. Normalizing the distributions by 
the number of GCs brighter than the turnover, we find that after 
correcting for incompleteness, px ~ 0.034. This is very similar to 
the directly observed value of px = 0.036 ± 0.05, so we conclude 
that any bias is very small. 

13 In both case, the uncertainties are computed assuming a 
Bernoulli distribution. 
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Fig. 12. — Color histograms for the full sample of globular clus- 
ters (open bars) and for the sample of globular clusters which coin- 
cide with an X-ray point source (filled bars). The curve is a normal 
kernel density estimate of the color distribution, which we denote 
by "S r (<?475 — ^850 )■ The tickmarks show the colors of individual 
GCs associated with X-ray point sources. 
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Fig. 13. — Histogram of (5475 — 2850) colors for 58 globular cluster 
in M87 which coincide with an X-ray point source. The curve shows 
a model color distribution of the form if> ~ lO0isir5-^ma)j/ j with 
(i = 0.77 (see text for details). The tickmarks show the colors of 
individual objects. 



To get a quantitative estimate of the dependence of px 
on metallicity, we model the color distribution of GCs 
which contain LMXBs as 

l/> ~ 10/3(9475-2=850) ^ (5) 

where "J/ is the color distribution for the full sample of 
GCs. To estimate ^, we use a kernel density estimate 
(Silverman 1986) which we denote by "i", and determine 
[3 via a maximum likelihood fit of the function 

ib ^ iQ^ 9475-2850 ^ (6) 

to the color distribution of the subsample of GCs that 
contain LMXBs. The result is (3 = 0.77 ± 0.22, and the 
distribution predicted by equation (6) is compared to the 
observed one in Figure IT31 To find the dependence of px 
on metallicity, it would be best to use an empirical deter- 
mination of the relation between (3475 — ^850) an d [Fe/H], 
but unfortunately no such relation is available in the lit- 
erature to the best of our knowledge. As a substitute, 
we use the models of Bruzual & Chariot (2003) to find 
the relation between (3475 — £850) and [Fe/H]. Using the 
data listed in Table |3 we find a best-fit linear relation of 
(#475 - z S5 o) ~ (0.38 ± 0.05) [Fe/H] + (1.62 ± 0.06), with 
an rms scatter of roughly 0.1 mag 14 . Using this relation, 
we find 

•0 ~ 10 (0.32±0.10)[Fc/H]^ 

~ (z/z )(°- 32±0 - 10 )* (7) 

14 The relation between [Fe/H] and (#475 — 2850 ) is slightly bet- 
ter described by a quadratic relation, (9475 — 2850) ~ (0.124 ± 
0.017)[Fe/H] 2 + (0.622±0.034)[Fe/H] + (1.620±0.015). If we use this 

relation, we would find that V ~ 10°- 0992 [ Fc / H ] 2 +°- 4976 [ Fc / H ] V. 
We prefer to use the linear relation because it adequately repre- 
sents the relation between [Fe/H] and (9475 — 2850 ) in the range 
where the GCs used to derive %j} lie (0.9 ^ (9475 — 2850) <^ 1-6) 
and because it allows us to cast ip simply in terms of a power law 
in Z. Although we do not attach any special physical significance 
to a power law form, it describes the observed trend with relative 
simplicity. 



We also reanalyzed the metallicity dependence of LMXBs 
in the GC system of M49, using the (V — I) colors pre- 
sented in Maccarone et al. (2002). Fitting a model of the 
form of equation (6), we find flvi = 1-2 ± 0.5. The em- 
pirical color-metallicity relation of Barmby et al. (2000) 
then gives 

il) VI ~ io( - 19±0 - 08 )[ F °/ H ^ y/ 

~ (z/z )(°- 19±ao8 )* v/ [ ' 

consistent with our findings for M87. The color distribu- 
tion for the M49 data, along with the best-fit model, is 
shown in Figure ITU 

5.3. Dependence on Luminosity 

In addition to metallicity, luminosity plays an impor- 
tant role in determining whether a given GC will contain 
a LMXB (in the sense that LMXBs reside preferentially 
in the most luminous clusters). This has been observed 
to be the case in NGC 1399 (Angelini et al. 2001), M49 
(Kundu et al. 2002), four early-type galaxies analyzed 
in Sarazin et al. (2003) and NGC 5128 (Minniti et al. 
2004). Given that our ACS observations of M87 define 
the deepest, most complete sample of GCs yet assembled 
for any galaxy, we now examine the dependence of px on 
luminosity in M87. 

In Figure IT31 we show the Zsscrband luminosity func- 
tion for the full sample of 1688 GCs within the ACS field 
of view (open histogram), along with the correspond- 
ing luminosity function for those 58 GCs which coin- 
cide with an X-ray point source (filled histogram). It 
is clear that, in agreement with previous findings, the 
LMXBs are associated preferentially with the brighter 
GCs. A two sample KS test rejects the hypothesis that 
the two datasets were drawn from the same parent dis- 
tribution and a Wilcoxon rank sum test rejects the hy- 
pothesis that the parent distributions share the same lo- 
cation; the respective p-values are pks — 7 x 10~ n and 

Pwilcox — 7 X 10 
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Fig. 14. — Histogram of (V— J) colors for globular clusters in M49 
that coincide with an X-ray point source, according to Maccarone 
et al. (2002). The curve shows a model color distribution of the 
form 4>vi ~ 10^v-/( v - J )^ VJ) with /3y7 = 1.2 (see text for details). 
The tickmarks show the colors of individual objects. 
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Fig. 16. — 2850 band luminosity function for globular cluster 
candidates with an associated X-ray point source. The solid and 
dashed curves are model distributions of the form if) ~ L a rf>, with 
a = 0.89 and a = 1, respectively. The tickmarks show the magni- 
tudes of individual objects. 
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Fig. 15. — zgso-band luminosity function for the full sample 
globular cluster candidates (empty bars) and for those candidates 
that coincide with an X-ray point source (filled bars) . The curve is 
a normal kernel density estimate of the luminosity function, which 
we denote by 4>(m). The tickmarks show the magnitudes of indi- 
vidual GCs associated with X-ray point sources. 



Both Kundu et al. (2002) and Sarazin et al. (2003) ar- 
gue that the data are consistent with the probability per 
unit luminosity being constant. To represent the lumi- 
nosity function of those GCs that coincide with LMXBs, 
4>(m), we will adopt a probability density of the form 

0(m) ~ L a t>(m) 

~ 10- 04Qm $(m) W 

where $(m) is the distribution of GC magnitudes. At 
this point, we could model $(m) parametrically (e.g., by 
using a Gaussian and taking care to model the effects 
of incompleteness). However, we prefer to approximate 



$(m) using a normal kernel density estimate (Silverman 
1986) since this approach does not require us to make any 
assumptions about the true parent distribution, while at 
the same time, the effects of incompleteness are taken 
into account. 

Denoting the kernel density estimate by $(m), we de- 
termine the parameter a via a maximum likelihood fit of 
the function 10~ OAam $(m) to the observed distribution, 
4>(m). The result is a = 0.89 ± 0.12, which is consistent 
with the results of Sarazin et al. (2003), who found that 
the luminosity function of GCs with associated LMXBs 
was consistent with the form (f>(m) ~ L$(m). In Fig- 
ure ^3 we show the luminosity function of GCs with 
an associated LMXB, along with the best-fit LMXB 
luminosity functions <fi ~ L 89 $(m) (solid curve) and 
<fi ~ LQ(m) (dashed curve). 

5.4. Dependence on Encounter Rates 

For the first time for a galaxy outside the Local Group, 
we are able to examine possible variations of px with 
GC structural parameters. The number of LMXBs is 
expected to depend on these parameters since the bina- 
ries responsible for the X-ray emission are likely to have 
a dynamical origin. That is to say, compact binaries in 
which one of the components is a neutron star are proba- 
bly not primordial in nature, but have likely formed as a 
result of the tidal capture of neutron star, or by exchange 
interactions with pre-existing binaries (Verbunt 2002). 

The encounter rates, T, for both tidal capture and ex- 
change interactions satisfy 

F oc ^ (10) 

v 

where po is the central mass density of the GC, r c is the 
core radius and v is the relative velocity of the encounter 
(Hut & Verbunt 1983). An estimate of v can be obtained 
from the velocity dispersion of the cluster, which in turn 
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is related to the central density and core radius via the 
virial theorem, a oc r cy /po. To be sure, other factors will 
affect the encounter rates for a given cluster, such as the 
rate at which binaries formed by these mechanisms are 
subsequently disrupted, the mass function in the cluster 
core, the binary fraction, and the period distribution of 
the binaries (Verbunt 2002). In this section, however, we 
concentrate on the explicit dependence of T on structural 
parameters; the role of other factors, which may them- 
selves depend on metallicity and structural parameters, 
will be examined in t |5.5l In the absence of more de- 
tailed information, we will first test the simple scenario 
in which px is proportional to 



(11) 



To compute L, we use the relations between core radius 
and half-light radius, loglZ(c) = log(rft/r c ), and between 
the dimensionless luminosity, log£(c) = log(L/j r 3 ), 
and King concentration parameter, c = log(r t /r c ), which 
are given in Appendix B of McLaughlin (2000). Here, jo 
is the central luminosity density, r^ is the half-light ra- 
dius, r c is the core radius, 15 r t is the tidal radius and L 
is the cluster's V-band luminosity. Explicitly, if m is the 
Z850-band magnitude, DM = 31.03 the distance modu- 
lus of M87 (Tonry et al. 2001), My,© = 4.84 the absolute 
V-band magnitude of the Sun (Vandenberg & Bell 1985), 
jo the central V-band luminosity density, and Yy the V- 
band mass-to-light ratio, then the central mass density, 
in units of M© pc -3 , is given by the expression 

^Q-OA(m-DM-M v ,e+c v ) 

Pa = Tvh = TCV ; — 775T37 (12) 

where the core radius is given by r c = rh/lZ and cy is the 
color term needed to convert from 2 850 -band to V-band 
luminosity. The V-band mass-to-light ratios of Galactic 
GCs are consistent with a constant value of Yy = 1.45 
in solar units (McLaughlin 2000), which we henceforth 
adopt for the M87 GCs. In general, the color term cy is 
ill-defined as it will depend on age and metallicity. We 
assume that the bulk of the GCs in M87 are old and 
coeval, as is suggested by spectroscopic and photomet- 
ric age measurements (Cohen, Blakeslee & Ryzhov 1998; 
Jordan et al. 2002; Kissler-Patig et al. 2002). Assuming 
a mean age of 13 Gyr (Cohen et al. 1998), we obtain 
the relation between ((7475 — 2:350) and (V — zsso) using 
the population synthesis models of Bruzual & Chariot 
(2003). For each (3475 — -zsso) color, we linearly interpo- 
late (or extrapolate, if necessary) to find the correspond- 
ing (V — zsso) at each of the metallicities in the Bruzual 
& Chariot (2003) models. The data used to define the 
color-color relation are listed in Table 

For the full sample of 1688 GCs, we find a mean en- 
counter rate of (logT) = 5.8 ± 0.02. For comparison, 
the mean encounter rate for the subsample of 58 GCs 
which coincide with a LMXB is (logT) = 6.73 ± 0.09. 
A two sample KS test rejects the hypothesis that the 
distribution of T parameters for these two samples arise 
from the same parent distribution and a Wilcoxon rank 
sum test rejects the hypothesis that they have the same 

5 Note that in the notation of McLaughlin (2000) the core radius 
is referred to as rg. 
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Fig. 17. — Distribution of (3475 — .zgso) color as a function of 
encounter rate, V = Pq' 5 t^, for the full sample of 1688 globular 
clusters (small crosses). Large circles indicate the 58 globular clus- 
ters which contain a LMXB. 



location; the respective p-values are pks = 3 x 10~ n 
and p W iicox — 3 x 5~ 15 . This constitutes the strongest 
evidence to date that encounter rates play a key role in 
determining px ■ 

This finding is also consistent with the results of Poo- 
ley et al. (2003) and Heinke et al. (2003). Pooley et al. 
(2003) find that T is the main factor in determining the 
number of close X-ray binaries in Galactic GCs. Specifi- 
cally, they find px oc T after restricting their analysis to 
the subset of GCs which contain bonafide LMXBs (see 
also Verbunt & Hut 1987), although their conclusions 
are hampered by the limited sample size. Heinke et al. 
(2003) find that the population of quiescent LMXBs in 
Galactic GCs is consistent with their dynamical origin 
as indicated by T. In M31, Bellazzini et al. (1995) 
have shown that the central density of GCs which host 
LMXBs is higher than the mean central density of M31 
GCs, which also points to the importance of encounter 
rates in determining the presence of LMXBs in GCs. 16 
Taken together, there seems to be little doubt that ex- 
change interactions and tidal captures in GCs are largely 
responsible for the production of LMXBs. 

In Figure El we plot the measured values of T against 
(ff475 — z 85o) color for the full sample of GCs (small 
symbols), along with the subsample of GCs which con- 
tain LMXBs (large symbols). There is a clear ten- 
dency for the latter GCs to have higher than average 
encounter rates, and it is apparent that no obvious cor- 
relation exists between T and GC color (recall from § 15.21 
that metallicity is an important factor in determining 
Px)- Since the encounter rate and metallicity are un- 
corrected, we can assume they are independent so that 
Pxocpi(I>a([Rs/H]). 

16 Note that the variations of V are driven mainly by po rather 
than r c - the former quantity varies by ~ 3 orders of magnitude and 
has a strong dependence on luminosity, while the latter quantity 
varies by ~ 1 order of magnitude and has a modest dependence on 
color or luminosity. 
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In M5.3I we showed that there is a strong correlation 
between px and GC luminosity. In Figure IThI we plot r 
versus zsso-band magnitude for the full sample pf GCs 
(small symbols), as well as for the those GCs which are 
associated with a LMXB (large symbols). There is a 
clear tendency for T to increase with increasing lumi- 
nosity, which raises the possibility that the trend dis- 
cussed in M5.3I is a consequence of a more fundamental 
correlation between T and luminosity. In fact, we can 
predict the luminosity distribution of GCs which contain 
LMXBs, </>(to), under the assumption that px oc T. Since 
there is no correlation between luminosity and color in 
the M87 GC system (Whitmore et al. 1995; Harris, Har- 
ris & McLaughlin 1998), we can safely ignore the metal- 
licity dependence of px when making this prediction. To 
find the behaviour of T as a function of zsso-band mag- 
nitude, m, we fit a cubic smoothing spline to the data 
in Figure ITB1 17 The resulting relation, r(m), is shown in 
FigureEDas the smooth curve. Using a kernel density es- 
timate, $(m), of the magnitude distribution for the full 
sample of GCs, the predicted distribution for the sub- 
sample of GCs that contain LMXBs should then satisfy 
<j){m) ~ r(m)$(m). 

This predicted distribution is shown as the dotted line 
in Figure EH We stress that no fitting has been done 
in making this comparison; the only assumption is that 
px oc r. The predicted distribution is in reasonable 
agreement with the observed distribution, although it 
seems to underpredict somewhat the number of faint 
GCs. A one-sample KS test gives pks = 0.07, which does 
not allow us to reject the hypothesis that the observed 
distribution is explained by the assumption px oc T with 
better than 99% confidence. 

Rather than assume px oc T, we now consider the de- 
pendence of px on a more general quantity, r Q , which 
satisfies 

r Q cx p%rl (13) 

where a is a free parameter that is determined by fit- 
ting to the observed magnitude distribution of those 
GCs which contain LMXBs. Such a form for the en- 
counter rate has been considered previously by Johnston, 
Kulkarni & Phinney (1992) 18 , who find a ~ 1.3 from an 
analysis of pulsars in Galactic GCs, and by Johnston & 
Vcrbunt (1996) who also find a ~ 1.3 from an analy- 
sis of low- luminosity X-ray sources in Galactic GCs. By 
adding this additional parameter, we are able to account 
for possible variations in other factors which may influ- 
ence the probability of forming LMXBs, such as varia- 
tions in the initial mass function (IMF) or systematic 
variations in the relative importance of tidal captures 
and binary-neutron star exchanges (c/. Grindlay 1996) 
. We would like r a to share with T the property of be- 
ing uncorrelated with (5475 — zsso) so as to be able to 

17 A smoothing spline minimizes over all functions / with con- 
tinous second derivatives a compromise between the fit and the 
smoothness of the form YKVi ~ f{ x i)) 2 + ^ J(f"( x )) 2 dx where 
{ x iy Vi} are the data and A controls the degree of smoothness and 
is chosen via cross-validation (Hastie & Tibshirani 1990; Green & 
Silverman 1994). 

18 Johnston et al. (1992) and Johnston & Verbunt (1996) assume 
r a oc Pq J M c oc Pn r c- Using the rough correlation r c oc Pq ' 2 
(McLaughlin 2000) it follows that a ~ ctj + 0.8 relates our a to 
theirs. 



CO - 




CO H 

19 20 21 22 23 24 25 26 
z 85 o [AB mag] 

Fig. 18. — Plot of zsso-band magnitude versus logT for all glob- 
ular clusters (small crosses). Large circles indicate those clusters 
that coincide with a LMXB. The curve shows a smoothing spline 
fit to the full dataset. 



consider them independent variables and thus to write 
Px = pi(r Q )p2([Fe/H]). We therefore define r Q as 

T a EE p « r 2 10 (0.46-a0.33)( S4 T5-* 8 5o) ) ( W ) 

where we have made use of the fact that r\ ~ 
1 Q(-o.46±o.06)( fl r 4 r5-2s5o) an( j pQ ^ 1 q(o.33±o.09)(s 4 75-.%5o), 

The latter expressions are best-fit relations obtained from 
our data. To determine a, we fit via maximum likelihood 
a function of the form (f> ~ r Q $(m) to the observed dis- 
tribution, where, as before, $ is a normal kernel density 
estimate of the GCs magnitude distribution and the be- 
haviour of po and as a function of magnitude has been 
obtained with smoothing splines. While formally it is 
more sound to include the color factor, we note that it 
has a negligible effect on the derived a since there is no 
color-luminosity correlation among M87 GCs. The best- 
fit value is a — 1.08 ± 0.12, the resulting cj) (the solid 
curve in Figure El shows very good agreement with the 
data . 

To summarize, we have shown that the encounter-rate 
based quantity T is an important factor in determining 
px- Moreover, we have shown that T, and the more 
general parameter T a , can account quantitatively for 
the observed magnitude distribution of GCs containing 
LMXBs. Whether one wishes to assign a more funda- 
mental role to r or r a , rather than luminosity, becomes 
to some extent a matter of taste. However, we believe 
that it is more correct to view px as being dependent on 
encounter-rate based quantities, as there are strong the- 
oretical arguments to support this interpretation. This 
is not true of the alternative view that px depends fun- 
damentally on cluster luminosity (i.e., more stars do not 
necessarily imply a more favorable environment for the 
production of compact binaries). 

5.5. Implications of the Derived Form for px 
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Fig. 19. — 2850-band luminosity function for globular clusters 
with an associated LMXB. The solid curve is a model distribution 
of the form <j> ~ r a= i.o8 < I > > while the dashed curve shows a model of 
the form cj> ~ The tickmarks show the magnitudes of individual 
objects. 



Combining the results of the previous sections, our es- 
timate for px is 

Px ~ T a=lm l0^^- Z8 ^ 

^ ^1.08^-1^00.87(5475-2850) (15) 

~ p\<l(Z/Z Q f^. 

Our goal in this section is to use this empirical relation to: 
(1) test the validity of the various physical mechanisms 
which have been proposed for the production of LMXBs 
in dense stellar environments; and (2) understand the 
origin of the observed metallicity dependence. 

There have been some theoretical suggestions as to why 
Galactic LMXBs might be more common in metal-rich 
GCs. For instance, Grindlay (1987) noted that, if the 
IMF depends on metallicity in such a way as to get flat- 
ter with increasing metallicity, then this would provide 
a larger population of massive stars (the progenitors of 
neutron stars in the LMXBs). A second, rather different, 
possibility is the suggestion by Bellazzini et al. (1995) 
that stars of higher metallicity will have larger radii and 
higher masses, thereby leading to enhancements in the 
tidal capture rates in metal-rich environments. Bellazz- 
ini et al. (1995) also note that such stars will more easily 
fill their Roche Lobes, further enhancing the number of 
LMXBs in metal-rich GCs. Recently, Maccarone et al. 
(2004) have suggested that irradiation induced winds on 
the donor star can explain the observed trend with metal- 
licity. We now examine these suggestions in detail. 

5.5.1. Dependence of the Number of Compact Stars per 
Unit Initial Mass on Metallicity 

Let us denote the number of neutron stars produced 
per unit initial mass by v, the fraction of these neutron 
stars that are retained by the cluster by /, and the total 
encounter rate (including both tidal captures and binary 
exchanges) by T t . Then 

r t cx fvpl- 6 T%{*2 + <?s)T (16) 



where po is the central mass density, r c is the core radius, 
T is the timescale over which the process lasts, and 02 
and 03 are the respective cross sections for tidal capture 
(two-body) and binary exchange (three-body) interac- 
tions (Johnston et al. 1992). Note that the cross-section 
for tidal capture scales with the radius, R, of the captur- 
ing star, while the binary exchange cross section scales 
with the semi-major axis, a, of the binary: i.e., o~2 cx R 
and o~ 3 tx a. 

Let us first consider possible metallicity-dependent 
variations in the fraction of neutron stars that are re- 
tained in GCs. Johnston et al. (1992) show that, for 
Galactic GCs, / varies by a factor of 3 <^ / <^ 5, with 
the precise value depending on the assumed velocity dis- 
tribution of newborn neutron stars. For the same sample 
of GCs, however, po varies by more than two orders of 
magnitude. Therefore, even though more massive GCs 
should certainly retain a greater fraction of their neutron 
stars, the net effect on T t is much smaller than that com- 
ing from the increase in the central density po, and we 
thus neglect it. The inclusion of this effect would steepen 
the dependence of T t on po. Note also that for GCs 
metallicity does not correlate with mass (e.g., McLaugh- 
lin & Pudritz 1996), so we expect / to be independent 
of metallicity. 

We follow usual practice in assuming that T is the 
same for the metal-rich and metal-poor subpopulations. 
Since the two subpopulations in M87 are observed to be 
roughly coeval (Cohen, Blakeslee & Ryzhov 1998; Jordan 
et al. 2002; Kissler-Patig et al. 2002), this assumption is 
reasonable, at least on average, although perhaps ques- 
tionable for individual clusters which may have suffered 
substantial dynamical evolution. 

Discounting a strong metallicity dependence in o~i or 
03 (see § I5.5.2|) , we are left with the result that the best 
candidate in T t to contain the metallicity dependence of 
px is v. Thus, we find that the observations point to 
a variation in the number of neutron stars formed per 
unit mass with GC metallicity, in the sense that more 
metal-rich GCs produce more neutron stars per unit ini- 
tial mass. Such a variation could be a consequence of 
variations in the IMF, as proposed by Grindlay (1987). 
But even for similar IMFs, metallicity will have impor- 
tant effect on stellar evolution which will affect the num- 
ber of neutron stars and black holes per unit mass. For 
instance, mass loss is thought to be greater for metal-rich 
stars and this will have a direct influence in the post- 
main-sequence evolution of high mass stars (Heger et al. 
2003). Thus, the metallicity dependence in v can plau- 
sibly arise through more than one process. Below we 
investigate on the possibility of IMF variations as the 
cause for the metallicity dependence, but it should be 
stressed that the result that v should be higher for GCs 
of higher metallicity is independent of the IMF being the 
cause of this dependence. 

In describing the IMF we will use a power law descrip- 
tion, with the number of stars with masses between m 
and m + dm given by N(m) cx mT x dm (where Salpeter 
is x — 2.35). This form is believed to be a good de- 
scription of the Galactic IMF only for stars with masses 
M J> 1M Q ; for lower masses a lognormal distribution 
appears to be a better description of the IMF (Chabrier 
2003). Nevertheless, we use the power-law form for 
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simplicity and to allow direct comparison with previous 
work. Our main focus below is the effect of IMF varia- 
tions in the relative number massive stars that end their 
evolution as neutron stars. Using a power law description 
this translates into a change in the IMF slope (assuming 
that all stars more massive than a certain value end up 
as neutron stars regardless of other factors) . 

If variations in v alone are responsible for the observed 
scaling of px with metallicity through IMF variations, 
then we can determine the behaviour that is required to 
produce the observed px- If all stars with m > m^s 
evolve to form neutron stars, then assuming that the 
minimum and maximum stellar masses are mi and m u , 
respectively, the number of neutron stars per unit mass, 
is, is given by (for x ^ 1, 2) 



v{x) 



1 



rn.„ 



NS 



(17) 



We parameterize the dependence of the IMF on metal- 
licity by assuming x depends linearly on [Fe/H], so 
dx/d[Fe/H] = A. We restrict the metallicity range to 
—2 < [Fe/H] < 0, which includes the vast majority of 
M87 GCs (Cohen et al. 1998). We assign x = -2.35 
to [Fe/H]= -1.0, assume m t = 0.08 Mq, m u = 100 M , 
m,vs = 8Mq, and then find A by minimizing the quan- 
tity 



(Cv(£,A) -£10°' 33C ) 2 dC, 



(18) 



where C and D are normalization constants. The re- 
sult of the minimization gives A = —0.3 ± 0.13 dex -1 , 
where the quoted uncertainty corresponds to the uncer- 
tainty in the metallicity dependence of px- Thus, the 
inferred metallicity dependence is fairly weak. This re- 
sult is in good agreement with the analysis of LMXBs 
in the Galactic and M31 GC systems presented by Bel- 
lazzini et al. (1995), who find that A ~ —0.4 dex -1 is 
necessary to be consistent with the ratio of LMXBs in 
metal-rich and metal-poor clusters. 

In their multivariate analysis of the Galactic GC sys- 
tem, Djorgovski et al. (1993) found that the intrinsic de- 
pendence of x on metallicity among Galactic GCs — af- 
ter removing the contribution of other important factors 
such as galactocentric distance, RgCi an d the distance 
from the Galactic plane, Zqp — is A w —0.5 dex -1 . 
Strictly speaking, these scaling relations are based on the 
present-day mass function in GCs, but it is reasonable to 
assume that this value of A reflects the dependence of the 
initial mass function on metallicity, as dynamical effects 
like cluster evaporation might be implicitly accounted 
for by the dependence on R gc and Z gp (Stiavelli et al. 
1991). 19 In any event, it is remarkable that our estimate 
for the required metallicity dependence of IMF slope is 
in good agreement with that found in the Galaxy. Thus, 
a difference in the IMF between the chemically distinct 
GC subpopulations remains a viable explanation for the 
observed metallicity dependence of px ■ 

19 The mass function in globular clusters will undoubtedly 
change due to dynamical effects, but as the stellar evolution 
time scales for M > 8Mq are shorter than dynamical evolution 
timescales, v should be sensitive mainly to the initial mass func- 
tion. 



5.5.2. Dependence of Stellar Radii on Metallicity 

Bcllazzini et al. (1995) argue that there is a second 
factor that can, in principle, increase px for metal rich 
GCs. According to these investigators, stars in high- 
metallicity GCs will have larger radii and masses than 
those in metal-poor GCs, leading to an enhancement in 
their tidal capture cross sections. To estimate the mag- 
nitude of the effect, Bellazzini et al. use the expression 
for the tidal capture rate, T TC , given in Lee & Ostriker 
(1986): 

T TC cx R 2 ~ T M T N. (19) 

Here R is the radius of the capturing star, M is its mass, 
N the total number of such stars, and the power-law in- 
dex, r = 1.07, is appropriate for our case (Lee & Ostriker 
1986). Using the stellar evolution models of Vanden- 
berg & Bell (1985), Bellazzini et al. (1995) find a ratio 
of capture rates for metal-rich and metal-poor stars of 
r\ = ^mrI^mp ~ 2.2. The enhancement is therefore 
comparable to the factor of three difference in px which 
we find for the two GC subpopulations in M87. 

However, the precise route to this ratio is not spelled 
out, and we were unable to reproduce their result. To 
estimate 77, we use isochrones of the Bergbush & Van- 
denberg (1992), adopt [Fe/H]= -0.47 and [Fe/H]=-1.82 
for the metal-rich and metal-poor GC subpopulations, 
respectively, and assume an age of 13 Gyr for both sub- 
populations. From equation 1191 we find rj ~ 2 if we 
compare two stars at the tip of the red giant branch, but 
this hardly constitutes a representative estimate for the 
whole population. The latter comparison might be mis- 
leading though, as the changes in the structure of a star 
as it ascends the red giant branch make the application 
of the cross section of Lee & Ostriker (1986) dubious. 
McMillan, Taam & McDermott (1990) find that the crit- 
ical impact parameter in units of the star radius is smaller 
for red giant stars than main sequence stars, and that it 
decreases as a star evolves through the red giant branch. 
In what follows we will neglect this effect, but note that 
inclusion of it would further reduce the magnitude of the 
effect advocated by Bellazini et al. (1995). 

A representative estimate of T TC may be obtained as 
follows. Assuming that the IMF is described by a power 
law, N(m) = m~ x dtn, we set 



where i runs over all tabulated masses. We thus weight 
the tidal capture rates by the expected number of stars 
at each mass. Using this approach, we find 77 ~ 1.1, with 
a modest dependence on the assumed value of x. Thus, 
given similar IMFs, this effect enhances the tidal capture 
probability of the metal-rich GCs by a negligible amount. 

A difference in IMF slope between the two GC sub- 
populations (see H5.5. 1[) has some effect on the estimated 
ratio. In this case, a correction factor — similar in form 
to the one in eauation ll7l with mjvs replaced by mi, and 
m u and m; set to the maximum and minimum masses 
in the isochrones — must be included to account for the 
fact that the number of stars per unit mass depends on 
x. Denoting such factor by n(x), the ratio is then 

n(xMR)rlf R (xMR) 



V 



n(xMp)Tjf P (xmp ) 



(21) 
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Taking x = 2.35 for the metal-poor population and 
x = 1.7 for the metal-rich (appropriate for the metallic- 
ity dependence of IMF slope in Galactic GCs, according 
to Djorgovski et al. 1993), we find rj ~ 1.3. If instead 
we set x = 1.35 for the metal-rich IMF, then the ratio 
increases to rj ~ 1.5. Thus, even under rather extreme 
assumptions for the metallicity dependence of IMF slope, 
the enhancement is insufficient to explain the observed 
factor of three difference in px- And, in any event, such 
a difference in x would result in a much larger enhance- 
ment in px through the increase in the relative numbers 
of neutron star progenitors (i.e., see H5.5.1[l . 

To summarize, we suggest that if the form of px is 
determined solely by dynamical processes, an increase 
in relative number of high-mass stars forming in metal- 
rich environments remains the only viable explanation 
for the observed metallicity dependence of px , although 
a full treatment of the problem should take into account 
the modest variations in tidal capture rates expected for 
stars of differing metallicity. 

5.5.3. Irradiation Induced Stellar Winds 

The discussion so far has assumed that Z has no ef- 
fect on intrinsic properties of LMXBs, such as their typ- 
ical lifetimes and luminosities, which might also produce 
the observed dependence of px on metallicity. Recently, 
Maccarone et al. (2004) proposed that irradiation in- 
duced stellar winds can explain the metallicity depen- 
dence of px ■ The basic mechanism is that radiation in- 
duced winds would be stronger in metal-poor donor stars 
due to less efficient metal line cooling and this would 
speed up the evolution of LMXBs in metal-poor clusters, 
leading to the observed trend assuming other processes 
such as the ones depicted above are not effective. They 
note that this mechanism may also explain the harder 
spectra observed in metal-poor Galactic LMXBs. 

Even though they do not provide a scaling relation that 
can be contrasted directly with the form we determine 
for px, they argue that the ratio of LMXBs for metal 
rich and metal-poor GCs (of metallicities Z r and Z p re- 
spectively) will scale roughly as (Z r /Z p )°' 3 ~ 0A . This ex- 
ponent in this scaling is very similar to the one we derive 
for px, and thus the observed form of px is certainly 
consistent with this scenario. 

5.5.4. Disruption and Hardening of Binaries? 

In terms of T, the dependence of px can be written 
px ~ r p^ n - 4 2 (Z/Z )°- 33 . It is worth reiterating the find- 
ings of q5.4l namely, that a ~ 1.1, which is reflected 
in the factor Pq 0A2 above. This result is in agreement 
with the findings of Pooley et al. (2003) find that the 
number of close X-ray binaries in Galactic GCs goes as 
N cx r°- 74±0 - 36 , which translates into a ~ 1.2 assum- 
ing a cx p^ 5 r c . Furthermore, the results are consistent 
with the findings of Johnston et al. (1992) and Johnston 
&: Verbunt (1996), which are not based on LMXBs, but 
rather on the statistics of pulsars and low-luminosity X- 
ray sources in Galactic GCs. All in all, these findings 
point to a scaling of px which is shallower than implied 
by the value of a = 1.5 in equation II II 

This reduction of a from its "expected" value of a — 
1.5, has potentially important implications for the forma- 
tion and evolution of LMXBs. For instance, one possible 



explanation for this weakening of the dependence of T t 
on po would be a change in the relative importance of 
o"2 and (T3 as a function of GC central density. Specif- 
ically, we can write the factor which includes the cross 
sections in T t (see equation I16|) as C = o~2 + C(Po)°~3- 
The factor C(i°o) could then account for the fact that, 
in denser clusters, stellar encounters would be more ef- 
fective in hardening wide binaries, thereby reducing 173. 
Since the stellar radii do not depend on GC structural 
parameters, we assume a 2. has no dependence on them. 
Then, from our best-fit form for px, we have C oc p$ ' , 
so that C(/°o) = (pPo° A — cr 2)/c r 3, where p is a constant. 
This scaling should serve as a constrain for predictions 
on the amount of binary hardening should this be the 
cause of the observed reduction on the expected value of 
a. 

Another possibility is that the destruction of binaries is 
responsible for the observed lessening of the dependence 
of r t on pq. Once formed, binaries will be destroyed at 
a rate Ar that satisfies (Verbunt 2002) 

ATtx/^V. (22) 
A simple model would then be that px is proportional 
to the ratio T/AT cx po r c ~ Po' 8r ci which is shallower 
but still marginally consistent with the observed, best-fit 
dependence of px on central density. 

6. SUMMARY AND CONCLUSIONS 

We have carried out the first detailed study of LMXBs 
in M87, using a catalog of 174 X-ray sources identified 
from deep Chandra/ACIS observations. All but ~ 20 
of these sources are expected to be LMXBs residing in 
M87. Combining the X-ray catalog with deep ACS imag- 
ing in the (74752:850 bandpasses for the central 11 arcmin 2 
of the galaxy, we have explored the connection between 
GCs and LMXBs. Our analysis is based on the largest 
sample of GC-LMXB associations currently available for 
any galaxy, and provides a first glimpse into the relation 
between compact accretors and their host GCs in M87. 

The luminosity function of X-ray sources is consis- 
tent with a single power law having an upper cutoff at 
Lx ~ 10 39 erg s _1 . Our reanalysis of data in the liter- 
ature (Kundu et al. 2002; Sarazin et al. 2001) reveals 
this also to be the case for M49 and NGC 4697; a similar 
conclusion was reached by Sivakoff et al. (2003) in their 
analysis of the LMXB luminosity functions in NGC 4365 
and NGC 4382. We conclude that, contrary to some pre- 
vious suggestions, there is no convincing evidence for a 
break in the luminosity function at Lx ~ 3 x 10 38 erg s _1 
(i.e., the Eddington luminosity corresponding to spher- 
ical accretion of ionized hydrogen onto the surface of 
a 1.4M Q neutron star). Given the sensitivity of this 
luminosity to the nature of the accretion process and 
the chemical composition of the accreted material, there 
seems to be no a priori reason to expect a sharp break 
in the luminosity function; indeed, we show through nu- 
merical simulations that the features identified by some 
previous researchers as breaks in the observed luminos- 
ity functions might be a consequence of the distribution 
of luminosities having an upper bound. These findings 
call into question the usefulness the luminosity function 
as a distance indicator, and cautions against using the 
inferred breaks to draw conclusions about black hole ac- 
cretors in early-type galaxies. If present, such black hole 
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accretors are better probed through studies of the spec- 
tral properties of the detected sources (e.g. Irwin et al. 
2003). The luminosity distribution of LMXBs could re- 
main a viable distance indicator if its form proves to be 
universal. The power law exponents we found by fitting 
truncated power laws are marginally consistent with a 
mean value of (7) = —1.78 ± 0.08 for the galaxies we 
considered; further studies of expanded samples should 
be able to test for variations in the luminosity function 
slope. 

In terms of LMXB formation efficiency in GCs, M87 
appears similar to other well-studied early-type galaxies. 
We find the percentage of GCs which contain LMXBs to 
be fx — 3.6 ± 0.5%, perfectly consistent with the values 
of 2 <^ fx <C, 4% found for a wide variety of early-type 
galaxies (Sarazin et al. 2003). The metal-rich GCs in 
M87 are observed to be 3±1 times more likely to contain 
LMXBs than the metal-poor GCs, consistent with pre- 
vious findings for M49 (Kundu et al. 2002). All in all, 
these results for LMXBs mirror other apparently "univer- 
sal" properties of GCs, most notably their near-Gaussian 
luminosity function (e.g., Harris 2001) and apparently 
constant formation efficiency (e.g., Blakeslee et al. 1997; 
McLaughlin 1999). Indeed, the constancy of fx, when 
coupled with the constant GC formation efficiency per 
unit baryon mass (Mclaughlin 1999), implies a constant 
LMXB formation efficiency per unit baryon mass in GCs. 
It would be interesting to investigate the bahaviour of the 
total number of LMXBs per unit baryon mass. This could 
have implications for the proposal that most LMXBs may 
form in GCs (White et al. 2002, Grindlay 1988), and that 
the presently observed populations of "field" LMXBs are 
the result of GC disruption and/or LMXB ejection via 
stellar encounters. As noted by White et al. (2002), 
support for this idea is provided by the observed scaling 
of the global LMXB X-ray luminosity to galactic optical 
luminosity, Lx, g iob/ L opt , with globular cluster specific 
frequency (van den Bergh & Harris 1981). If all LMXBs 
are formed in GCs, then a constant LMXB formation ef- 
ficiency would arise naturally; on the other hand, if there 
are separate populations of field and GC LMXBs, with 
different origins, it becomes harder to explain, as the 
LMXB population arising from field stars would have to 
know about the fraction of baryons that are not in the 
form of stars. In M87, the observed similarity between 
the luminosity distributions of field and GC LMXBs is 
broadly consistent with a scenario in which most LMXBs 
form in GCs. 

In agreement with previous findings based on smaller 
samples of LMXB-GC associations, we find that both 
GC metallicity and luminosity are important factors in 
determining the presence of LMXBs in GCs (Kundu et al. 
2002). Furthermore, we have been able to demonstrate 
that the probability, px, that a given GC will contain a 
LMXB depends sensitively on the parameter T = Po 5 r 2 , 
which is proportional to the tidal capture and binary- 
neutron star exchange rates within the host GC. This 
constitutes the strongest evidence to date that these dy- 
namical processes are responsible for the formation of the 
bulk of LMXBs in GCs. 

Working from the subsample of GCs which contain 
LMXBs and which have colors and magnitudes measured 
from our deep ACS images, we have explored the scaling 



of px with a variety of GC structural and photometric 
parameters. We confirm the previously identified depen- 
dence of px on GC metallicity and luminosity, but ar- 
gue that the observed luminosity dependence arises as 
a result of the enhanced encounter rates for more lumi- 
nous clusters (mainly because the central mass densities 
of GCs increase with increasing luminosity; McLaughlin 
2000). Considering the dependence on structural param- 
eters more fundamental, our preferred expression for px 
is then 

px oc Tp^ A2±0 - 11 (Z/Z G) f^ ±a -\ (23) 

The metallicity dependence in this scaling relation - 
which translates into the aforementioned factor of three 
enhancement in px for metal-rich GCs relative to their 
metal-poor counterparts — has in the past been pro- 
posed to be a result of metallicity-dependent variations 
in the IMF (Grindlay 1987), a consequence of irradia- 
tion induced stellar winds (Maccarone et al. 2004) , or of 
an enhancement in tidal-capture rates due to the larger 
radii of metal-rich stars (Bellazzini et al. 1995). 

We critically examine the viability of these mechanisms 
in light of the new observation constraints for M87, and 
find that previous studies have likely overestimated the 
importance of the latter mechanism. Assuming a uni- 
versal power law IMF, our calculations suggest a typical 
enhancement of ~ 10% for the metal-rich GCs due to 
this radius-metallicity dependence, far smaller than the 
observed factor-of-three difference. Only by allowing the 
IMF to vary between the chemically distinct GC sub- 
populations it is possible to produce such enhancements, 
but even in this case, the increase in px is driven mainly 
by the increased number of neutron star progenitors in 
metal-rich environments. On the other hand, the de- 
pendence of IMF slope x on metallicity which is needed 
to account for the observed metallicity dependence of 
px is found to be dx/d[Fe/R] = -0.3 ± 0.13 dex" 1 . A 
variation of the IMF slope with metallicity produces the 
metallicity dependence in px by increasing the number 
of compact stars per unit initial mass for metal-rich GCs. 
The need for an increased number of compact stars for 
higher Z is independent of the particular form of the 
IMF. We conclude that the only viable dynamical means 
of accounting for the observed metallicity dependence of 
px appears to be a relative enhancement in the number 
of neutron stars in metal-rich GCs. It is possible that 
intrinsic properties of LMXBs, unrelated to dynamical 
properties of the host GC, are affected by Z and these 
in term can affect the form of px- One such mecha- 
nism, irradiation induced winds, has been recently pro- 
posed (Maccarone et al. 2004) and it is consistent with 
the observed form of px- Further detailed studies of 
the proposed mechanism and contrasting its predictions 
with observations should shed light on which mechanism 
- dynamical or intrinsic - determines the form of px ■ 

What emerges from our observations is that a sim- 
ple dynamical picture — namely, the capture of neu- 
tron stars by single or binary stars within GCs, as fully 
expressed in equation (|lfi|) — can account quantitatively 
for the observed scaling of px with structural parameters 
and metallicity. A detailed investigation of the observed 
behaviour of px should be now undertaken, along with a 
comparison to the results of numerical simulations which 
probe the formation of compact binaries in dense stellar 
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environments and incorporate realistic stellar structure 
and evolution models. It would be particu larly useful to 
examine the possibility, discussed in H5.5.41 that the effec- 
tive encounter rates are reduced by a factor Ar ~ Pq° A - 
Thus, our findings for the LMXB population in M87 may 
be evidence for the ongoing disruption of binary stars in 
dense environments, or a reduction in the binary-neutron 
star exchange rates due to the hardening of compact bi- 
naries via close encounters (e.g., Hut et al. 1992). 

To date, studies of the connection between LMXBs and 
GCs based on Chandra observations of external galaxies 
have been hindered by the lack of high-quality optical 
data needed to characterize their GC systems. The ACS 
Virgo Cluster Survey (Cote et al. 2004) will offer a nearly 
complete census of GCs within the central regions of one 
hundred early-type galaxies in the Virgo Cluster. The 
measurement of metallicities, luminosities, perhaps most 
importantly, structural parameters, for the many thou- 
sands of GCs which will be detected in the course of this 
survey will allow refinements to the detailed form of px , 
including an exploration of its possible dependence on 
galaxy environment. 
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X-ray Point-sources in M87 
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-0.09 


-0.84 


no 




2 


cxou 


J123034.84 


-122527 


2.98 


4.29 


0.17 


-0.13 


no 




3 


cxou 


J123034.84 


-122215 


12.67 


18.44 


0.68 


0.67 


no 




4 


cxou 


J123035.84 


-122341 


3.11 


4.42 


0.07 


-0.37 


no 




5 


cxou 


J123036.04 


-122433 


5.41 


7.71 


0.26 


0.11 


no 




6 


cxou 


J123036.24 


-122532 


13.29 


19.11 


0.19 


-0.04 


no 




7 


cxou 


J123036.64 


-122213 


3.77 


5.45 


-0.15 


-0.03 


no 




8 


cxou 


J123036.84 


-122459 


2.78 


3.96 


0.15 


0.48 


no 




9 


cxou 


J123038.94-122304 


5.31 


7.49 


0.00 


-0.52 


no 




10 


cxou 


J123039.94-122550 


2.57 


4.21 


-0.12 


-0.31 


no 




11 


cxou 


J123040.34-122509 


6.37 


8.97 


-0.17 


-0.36 


no 




12 


cxou 


J123040.9+122320 


7.89 


11.14 


0.17 


-0.14 


no 




13 


cxou 


J123041. 04-122403 


18.34 


25.21 


0.07 


-0.17 


no 




14 


cxou 


J123041. 2+122450 


4.14 


5.81 


-0.29 


0.09 


no 




15 


cxou 


J123041. 24-122308 


1.77 


2.51 


-0.45 


-1.00 


no 




16 


cxou 


J123041. 64-122203 


3.48 


4.81 


-0.06 


0.22 


no 




17 


cxou 


J123041.64 


-122600 


10.45 


14.71 


-0.02 


-0.25 


no 




18 


cxou 


J123041.74 


-122439 


18.44 


25.78 


0.09 


-0.09 


no 




19 


cxou 


J123042.04 


-122626 


1.70 


2.35 


-0.52 


-0.85 


no 




20 


cxou 


J123042.04 


-122450 


2.35 


3.29 


-0.35 


-0.65 


no 




21 


cxou 


J123042.34 


-122156 


6.66 


9.23 


0.03 


-0.05 


no 




22 


cxou 


J123042.74 


-122135 


4.34 


8.62 


-0.04 


-0.64 


no 




23 


cxou 


J123043.04 


-122525 


3.36 


4.71 


0.17 


0.07 


no 




24 


cxou 


J123043.14 


-122502 


7.31 


10.20 


-0.13 


-0.04 


no 




25 


cxou 


J123043.44 


-122422 


5.78 


8.01 


0.18 


-0.00 


no 




26 


cxou 


J123043.44 


-122746 


2.85 


4.06 


0.28 


-0.34 


no 




27 


cxou 


J123043.54 


-122346 


9.84 


14.35 


0.12 


-0.10 


no 




28 


cxou 


J123043.74-122429 


3.77 


5.24 


0.68 


0.25 


no 




29 


cxou 


J123043.74-122418 


3.26 


4.51 


-0.16 


-0.60 


no 




30 


cxou 


J123044.04-122307 


1.97 


2.74 


0.47 


0.21 


yes 


no 


31 


cxou 


J123044.1+122456 


5.70 


7.88 


0.06 


-0.22 


no 




32 


cxou 


J123044.24-122312 


2.82 


3.90 


-0.26 


-0.40 


yes 


globular 


33 


cxou 


J123044.2+122134 


14.09 


19.24 


-0.17 


-0.22 


no 




34 


cxou 


J123044.2+122209 


20.72 


28.49 


0.09 


-0.06 


yes 


no 


35 


cxou 


J123044.54 


-122254 


3.94 


5.47 


0.62 


0.05 


yes 


no 


36 


cxou 


J123044.54 


-122450 


3.97 


5.47 


0.65 


0.02 


yes 


globular 


37 


cxou 


J123044.64 


-122140 


27.39 


38.06 


0.17 


-0.25 


no 




38 


cxou 


J123044.64 


-122201 


23.83 


32.65 


0.27 


-0.06 


yes 


globular 


39 


cxou 


J123044.64 


-122328 


5.20 


7.39 


0.28 


-0.10 


yes 


globular 


40 


cxou 


J123044.74 


-122434 


49.66 


68.86 


-0.20 


-0.55 


yes 


non globular 


41 


cxou 


J123044.94 


-122404 


8.81 


12.26 


0.44 


0.10 


yes 


no 


42 


cxou 


J123044.94 


-122436 


3.24 


4.49 


-0.44 


-0.38 


ves 


globular 


43 


cxou 


J123045.04 


-122317 


5.03 


6.77 


0.38 


0.20 


yes 


no 


44 


cxou 


J123045.24-122425 


5.36 


7.48 


0.48 


0.45 


yes 


globular 


45 


cxou 


J123045.34-122352 


1.98 


2.75 


0.60 


0.55 


yes 


globular 


46 


cxou 


J123045.44-122702 


5.33 


7.33 


0.02 


0.13 


no 




47 


cxou 


J123045.44-122519 


3.82 


5.54 


0.42 


-0.10 


no 




48 


cxou 


J123045.4+122329 


2.93 


4.19 


0.21 


-0.65 


yes 


no 


49 


cxou 


J123045.54-122412 


3.34 


4.66 


0.63 


0.31 


yes 


globular 


50 


cxou 


J123045.54-122450 


4.44 


6.11 


-0.17 


-0.86 


yes 


no 


51 


cxou 


J123045.74-122409 


6.95 


9.71 


0.30 


-0.02 


yes 


globular 


52 


cxou 


J123045.84 


-122134 


3.64 


5.07 


0.38 


-0.10 


no 




53 


cxou 


J123045.84 


-122336 


3.36 


5.35 


-0.12 


-0.30 


yes 


no 


54 


cxou 


J123045.94 


-122408 


4.42 


6.17 


0.29 


-0.12 


yes 


no 


55 


cxou 


J123045.94 


-122125 


4.96 


7.68 


-0.15 


-0.12 


no 




56 


cxou 


J123046.24 


-122328 


18.40 


26.27 


-0.03 


-0.34 


yes 


globular 


57 


cxou 


J123046.24 


-122349 


4.58 


6.36 


-0.03 


0.17 


yes 


globular 


58 


cxou 


J123046.34 


-122207 


3.00 


4.12 


0.63 


-0.62 


yes 


no 


59 


cxou 


J123046.34 


-122323 


14.49 


20.62 


0.17 


-0.01 


yes 


globular 


60 


cxou 


J123046.34 


-122432 


3.96 


5.50 


-0.80 


-0.13 


yes 


globular 


61 


cxou 


J123046.34 


-122441 


3.47 


4.79 


-0.07 


-1.00 


yes 


globular 


62 


cxou 


J123046.54 


-122450 


6.75 


9.33 


-0.17 


-0.16 


yes 


globular 


63 


cxou 


J123046.64 


-122223 


5.42 


7.48 


-0.25 


-0.37 


yes 


globular 


64 


cxou 


J123046.74-122402 


9.00 


12.57 


0.15 


-0.22 


yes 


see note 1 


65 


cxou 


J123046.74-122237 


3.41 


4.71 


-0.51 


-0.24 


yes 


globular 


66 


cxou 


J123046.74-122453 


4.84 


6.70 


0.00 


-0.12 


yes 


globular 


67 


cxou 


J123046.74-122150 


4.98 


6.82 


-0.09 


-0.51 


yes 


globular 


68 


cxou 


J123046.8+122305 


8.20 


11.19 


-0.34 


-0.04 


yes 


no 


69 


cxou 


J123046.94-122259 


3.74 


5.08 


0.05 


-0.15 


yes 


no 


70 


cxou 


J123047.04-122518 


2.27 


4.14 


-0.28 


-0.25 


no 




71 


cxou 


J123047.04 


-122500 


4.59 


6.37 


-0.19 


-0.52 


yes 


no 


72 


cxou 


J123047.14 


-122415 


99.35 


139.27 


-0.07 


-0.65 


yes 


see note 2 


73 


cxou 


J123047.34 


-122308 


14.80 


20.10 


0.08 


-0.23 


yes 


no 


74 


cxou 


J123047.54-122126 


9.66 


13.54 


-0.04 


-0.10 


no 





JORDAN ET AL. 

TABLE 1 — Continued 



ID 




IAU name 




count rate 
10~ 4 cts s -1 


(10 J7 erg s L ) 


H21 b 


H31 b 


ACS 


opt. counterpart? 


7^ 


fYOTT 


J123047.54 


-122621 


9 QQ 
z.yy 


A Ofi 


9fi 
-U.ZU 


1 fi 


no 


' ' ' . 


7fi 
t 


PYOT T 


J123047.54 


-122423 


7 80 
/ .OU 


1 80 

iu.oy 


14 
-U. 14: 


A(\ 

-U.4:U 


yes 


globular 


77 


fYOTT 

L AU VJ 


J123047.64 


-122351 


Q 1 8 

y . -Lo 


1 9 80 
-LZ.OU 


n i n 

U. 1U 


fiO 

-u.uu 


yes 


no 


78 
/ 


PYfiT T 

L.AU VJ 


J123047.64 


-122234 


A ^A 


« 9^ 


fl 1 ^ 
U. 10 


04 

-U.U4: 


yes 


no 


7Q 


fYOTT 

L AU U 


J123047.64 


-122220 


7 90 


Q 09 

y.yo 


09 

-u.uo 


11 
-U. 1 1 


yes 


globular 




PYfiT T 

UAU U 


J123047.74 


-122334 


3^ D7 


48 91 


14 

U. 14: 


09 
-U.UO 


yes 


globular 




PYflT T 
LAU VJ 


J123047.84 


-122404 


7 ^o 
1 .ou 


10^1 
1U.01 


n 1 7 

U. 1 i 


90 
-U.ZU 


yes 


globular 


oz 


fYOTT 

L AU U 


J123047.84 


-122401 


A 88 


8^ 


-0 12 


47 

-U.4: / 


yes 


globular 


89 


fYOTT 
L AU U 


J123047.9+122618 


3 ^A 


4 8ft 


n 94 

U.Z4: 


10 
-U. 1U 


no 




84 
04: 


PYOT T 
LAU U 


J123048.04-122431 


1 71 
O. 1 1 


^14 


94 

-U.Z4: 


^9 
-U.OO 


yes 


' ' ' , 
globular 


8^ 


PYOTT 


J123048.3+122455 


4 54 


99 


-0 14 


1^ 
U. IO 


yes 


globular 


8fi 
OO 


v. au u 


J123048.3+122415 


O.Ol 


7 79 


fl 9^ 
U.00 


07 
-U.U / 


yes 


globular 


87 


PYOTT 

vAU U 


J123048.34-122438 


4 74 


^7 


fl 07 
U.U i 


01 
-U.U1 


yes 


globular 


88 


UAU U 


J123048.74-122620 


A QQ 

^.yy 


7^ 


fl ^7 
U.O ( 


99 
-U.OZ 


no 




8Q 


PYPiTI 
UAU U 


J123048.74-122517 


3 £Q 

o.oy 


^ 08 

o.yo 


Ofi 
-U.UU 


91 
-U.Ol 


no 


' ' ' , 


on 
yu 


PYflT T 
UAU VJ 


J123048.74-122414 


Q ^9 

y.oz 


1 9 91 
10. Ol 


fl 91 
U.Ol 


14 

U. 14: 


yes 


globular 


y i 


PYPiT T 

UAU U 


J123048.84 


-122347 


1Q.O0 


99 81 
ZZ.ol 


07 
U.U 1 


48 

-U.4:0 


yes 


globular 


Q9 

yz 


PYOTT 

UAU VJ 


J123048.84 


-122313 


i "\ 9n 


1 8 fifi 


17 
U. 1 1 


ZL8 

~U.4:0 


yes 


no 


09 

yo 


PYPiTT 
UAU VJ 


J123048.94 


-122343 


1 £. RQ 
ID.Dy 


99 97 


00 

-u.uy 


74 

~U. 1 4: 


yes 


no 


OA 

y^i 


PYOT T 
UAU U 


J123049.04 


-122405 


Q fil 

y.oi 


1 9 48 
10. 


n 17 

U. 1 ( 


^9 
-U.OO 


yes 


globular 


yo 


PYOTT 

UAU VJ 


J123049.14 


-122159 


Q 80 
y.ou 


1 9 ^9 


91 
-U.Ol 


-0 12 


yes 


n .° , . 


OR 

yo 


PYflT T 

UAU U 


J123049.14 


-122445 


9 7fi 


9 89 
o.oz 


n 48 

U.4:0 


40 
U.4:y 


yes 


globular 


Q7 


PYOTT 

UAU VJ 


J123049.14 


-122604 


Q3 °,7 


1 9fi 71 
1ZU. I 1 


04 

-U.U4: 


^0 

-u.ou 


no 


',' ' . 


Q8 
yo 


PYOTT 

UAU U 


J123049.14 


-122308 


°, 71 


^ 87 


n 99 

u.oz 


-0 12 


yes 


glo bulcir 


yy 


PYOTT 

UAU U 


J123049.24 


-122334 


7Q Ofi 
t y.uu 


1 OQ ZL7 


-0 22 


80 

-u.ou 


yes 


globul&r 




PYOTT 

UAU VJ 


J123049.54 


-122355 


12 17 


17 05 


07 
U.U i 


^ 
-u.oo 


yes 


globuleir 


1 01 

1U1 


PYflT T 

UAU U 


J123049.64 


-122333 


9£ A 9 


9fi (\0 
OD.DZ 


11 
U. 1 1 


99 
-U.OO 


yes 


globular 


102 


PYOTT 

UAU U 


J123049.64 


-122353 


3 fi7 


^14 

O- 14: 


U.1U 


-0 42 


yes 


globulcir 


1 09 


PYOTT 

UAU VJ 


J123049.74-122351 


i ^ i o 

It). 1U 


21 12 


10 

u. iy 


90 

-u.ou 


yes 


globular 


1 04 
1U4: 


UAU VJ 


J123049.84-122402 


99 R9 


91 70 


99 

u.zz 


99 

-u.oz 


yes 


globular 


i o^ 

1UO 


UAU VJ 


J123049.84-122216 


9 38 


9 97 
O.Z ( 


-0 21 


fiQ 

-u.uy 


yes 


no 




fYOTT 
UAU VJ 


J123049.84T22436 


8 AQ 


1 1 79 
11.(0 


08 
-U.Uo 


90 
-U.ZU 


yes 


no 


1 07 

1U f 


fYOTT 

UAU VJ 


J123049.94-122740 


A 37 


u.oy 


n 97 

U.Z i 


09 
u.uz 


no 




1 08 


fYOTT 

UAU U 


J123050.0+122400 


1 7 fi7 


9/L 77 




97 

-U.Z / 


yes 


I K 1 

globular 


1 00 

iuy 


fYOTT 
UAU U 


J123050.14-122251 


i n Q8 
lu.yo 


1 4 07 
14:. y / 


04 

U.U4: 


4^ 

-U.4:0 


yes 


globular 


110 


fYOTT 

UAU VJ 


J123050.14 


-122301 


93 81 


99 ^0 

oo. ou 


10 
U. 1U 


9^ 

-u.oo 


yes 


globular 


ill 
i ± i 


fYOTT 
UAU VJ 


J123050.24 


-122608 




« 97 
U.Z ( 


11 
-U. 1 1 


47 
-U.4: / 


no 




112 


fYOTT 

UAU VJ 


J123050.34 


-122128 


3 39 


4 


n ^8 
u.oo 


01 
U.U1 


no 




no 


fYOTT 

UAU VJ 


J123050.34 


-122332 


5 65 


7 84 


n 99 

u.zo 


94 

-U.04: 


yes 


I K 1 

globular 


114 
114: 


fYOTT 
UAU VJ 


J123050.44 


-122212 


^ 70 
o. / u 


7 89 
/ .OZ 


99 

-u.oz 


94 

-U.04: 


yes 


globular 


1 1 ^ 

llO 


fYOTT 

UAU VJ 


J123050.54 


-122356 


1 8 38 
10.00 


9^ 7^ 
ZO. ( o 


17 

U. 1 1 


A(\ 

~U.4:U 


yes 


n .° , . 


1 1 (\ 
HO 


fYOTT 
UAU U 


J123050.54 


-122435 


3 11 


4 90 

4:.zy 


70 

-u. i y 


flA 

-U.04: 


yes 


globular 


117 
1 1 f 


fYOTT 

UAU VJ 


J123050.84 


-122502 


A3 1 ^ 


91 

uy .01 


12 


98 
-u.zo 


yes 


no 


118 

no 


fYOTT 

UAU VJ 


J123050.84 


-122411 


AQ 


04 

y.U4: 


99 
-u.oz 


11 
-U. 1 1 


yes 


n .° , . 


1 1 Q 

i iy 


fYOTT 
UAU U 


J123051. 1+122242 


3 ^9 


4 77 

4:. I ( 


n 4i 


00 

-u.uy 


yes 


globular 


120 


fYOTT 

UAU U 


J123051. 84-122159 


^ 47 


7 ZLQ 
( .4:y 


U.U i 


22 


yes 


no 


191 


fYOTT 
UAU U 


J123051. 84-122247 


9 8Q 

z.oy 


9 04 

o.y4: 


U. 10 


99 

u.zo 


yes 


no 


122 


fYOTT 

UAU U 


J123051. 84-122911 


1 41 


9 19 

z. io 


-0 12 


1 00 
-l.UU 


no 




123 


fYOTT 

UAU VJ 


J123052.54-122533 


A QO 


7^ 

u. / o 


n ^fi 
u.oo 


24 


no 




1 94 
1Z4: 


fYOTT 
UAU VJ 


J123052.64-122323 


3 Q^ 

o.yo 


^ 40 

o.4:y 


1 00 
l.UU 


1 00 
l.UU 


yes 


no 


125 


fYOTT 

UAU U 


J123052.74-122336 


7 1 Q 


1 flfi 
1U.UU 


n 98 
u.oo 


97 
U.Z i 


yes 


globular 


1 OR 
1ZO 


fYOTT 
UAU U 


J123052.94 


-122547 


7 3A 


00 

y.yy 


04 

-U.U4: 


91 
-U.Z1 


no 




1 97 


fYOTT 
UAU U 


J123053.04 


-122244 


3 Q7 

o.y / 


^ 4fl 


17 
-U. 1 / 


88 

-u.oo 


yes 


no 


128 


fYOTT 

UAU U 


J123053.04 


-122535 


4 27 


^ 84 

0.04: 


n 9^ 
u.zo 


98 
-u.zo 


no 




1 9Q 

izy 


fYOTT 
UAU U 


J123053.04 


-122208 


A AP, 


« 07 
U.U i 


fi9 
-U.OO 


18 
-U. IO 


yes 


no 


1 90 


fYOTT 

UAU VJ 


J123053.24 


-122356 


28 14 


90 40 

oy .4:U 


10 
U. 1U 


90 
-U.ZU 


yes 


globular 


191 
101 


fYOTT 
UAU VJ 


J123053.44 


-122556 


3 i fi 

O. ±0 


4 9£ 

4:.ZU 


n 99 
u.zz 


98 
U.ZO 


no 




1 99 


fYOTT 
UAU U 


J123053.64 


-122237 


A 01 


^ 44 

0.4:4: 


49 

-U.4:Z 


07 
-U.U / 


yes 


non globular 


1 99 
-LOO 


fYOTT 

UAU VJ 


J123053.74 


-122448 


2 91 


9 00 

o.yy 


n 9^ 

U.ZO 


^il 
-U.Ol 


yes 


n .° , . 


1 94 
104: 


fYOTT 
UAU U 


J123053.94 


-122430 


A 73 


fl ^1 
U.01 


49 

-U.4:Z 


19 
-U. 1Z 


yes 


globular 


1 9^ 

-LOO 


fYOTT 

UAU VJ 


J123053.94 


-122544 


^Q 
u.oy 


8 Qfi 

o.yu 


n 9^ 

U.ZO 


-0 44 


no 


' ' ' . 


1 9fi 
-LOU 


fYOTT 

UAU U 


J123054.44 


-122302 


7 01 


1 07 
1U.U I 


n 09 
u.uz 


14 


yes 


globular 


1 97 
-LO ( 


fYOTT 

UAU VJ 


J123054.64-122218 


1 Q7 
i .y / 


9 R7 

Z.U ( 


99 
-u.oz 


17 

U. 1 i 


yes 


no 


1 98 
-LOO 


fYOTT 

UAU U 


J123054.74-122222 




10 12 


n 98 

u.zo 


-0 22 


yes 


globular 


1 9Q 

ioy 


fYOTT 
UAU VJ 


J123054.84-122231 


3 ZL9 


4 fi^ 

4:. DO 


n 99 
u.zo 


A(\ 

-U.4:U 


yes 


no 


140 


cxou 


J123054.94-122247 


4.22 


5.98 


-0.04 


-0.46 


yes 


no 


141 


cxou 


J123054.94-122538 


22.17 


30.24 


0.07 


-0.12 


no 




142 


cxou 


J123054.94-122438 


8.14 


11.20 


0.27 


0.16 


yes 


globular 


143 


cxou 


J123055.04-122504 


2.83 


4.05 


0.57 


-0.76 


yes 


no 


144 


cxou 


J123055.24-122340 


11.83 


16.63 


-0.03 


-0.03 


yes 


no 


145 


cxou 


J123055.3+122256 


3.41 


5.23 


-1.00 


-0.50 


yes 


globular 


146 


cxou 


J123055.44 


-122314 


3.81 


5.31 


-0.34 


0.01 


yes 


globular 


147 


cxou 


J123055.44 


-122342 


10.93 


15.35 


0.15 


0.06 


yes 


globular 


148 


cxou 


J123055.54 


-122615 


2.64 


3.56 


0.24 


-0.13 


no 




149 


cxou 


J123056.2+122526 


2.49 


3.40 


1.00 


1.00 


no 
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TABLE 1 


— Continued 










ID 




IAU name 


COUIlt TcLlC 


' J .X 

('\ fl 37 pro- ^\ 
\ lu erg S ) 


H21 b 


H31 b 


ACS 


opt. counterpart? 


150 


cxou 


J123056.2+122447 


8.26 


11.79 


0.12 


-0.02 


yes 


globular 


151 


cxou 


J123056.3+122211 


11.71 


15.95 


0.12 


-0.18 


yes 


globular 


152 


cxou 


J123056.4+122448 


4.65 


7.39 


0.02 


-0.26 


yes 


globular 


153 


cxou 


J123056.7+122259 


6.44 


8.93 


0.06 


-0.50 


yes 


no 


154 


cxou 


J123057.2+122122 


14.09 


19.44 


0.13 


-0.04 


no 




155 


cxou 


J123057.9+122220 


3.92 


5.37 


0.48 


0.34 


no 




156 


cxou 


J123058.0+122844 


1.57 


2.31 


0.24 


0.13 


no 




157 


cxou 


J123058.1+122104 


9.51 


13.16 


0.20 


-0.31 


no 




158 


cxou 


J123058.5+122222 


10.03 


13.77 


0.13 


-0.25 


no 




159 


cxou 


J123059.0+122259 


3.47 


4.83 


0.13 


-0.59 


no 




160 


cxou 


J123059.5+122038 


38.97 


54.50 


-0.23 


-0.51 


no 




161 


cxou 


J123059.6+122509 


25.13 


34.53 


0.02 


-0.36 


no 




162 


cxou 


J123100.1+122232 


4.58 


6.58 


-0.16 


-0.51 


no 




163 


cxou 


J123100.2+122138 


3.37 


4.69 


-0.00 


-0.45 


no 




164 


cxou 


J123100.3+122417 


14.62 


20.30 


0.27 


0.08 


no 




165 


cxou 


J123100.6+122021 


1.30 


3.12 


-0.42 


-0.37 


no 




166 


cxou 


J123102.4+122040 


1.52 


2.14 


0.36 


0.15 


no 




167 


cxou 


J123102.6+122121 


2.21 


3.12 


0.18 


-0.70 


no 




168 


cxou 


J123102.6+122411 


5.06 


7.05 


0.04 


-0.52 


no 




169 


cxou 


J123103.3+122107 


3.71 


5.24 


-0.38 


-0.21 


no 




170 


cxou 


J123103.5+122305 


2.30 


3.22 


-0.17 


-1.00 


no 




171 


cxou 


J123104.1+122335 


1.71 


2.39 


-0.20 


-1.00 


no 




172 


cxou 


J123104.5+122156 


2.97 


4.80 


0.21 


0.14 


no 




173 


cxou 


J123105.5+122731 


2.35 


3.38 


0.28 


-0.30 


no 




174 


cxou 


J123108.9+122618 


1.03 


1.53 


0.27 


0.27 


no 





a Lx values are in the 0.3— 10 keV energy band and are obtained assuming a power law spectral shape with power law index 
k — 1.64 and a Galactic hydrogen column density Nh = 2.5 x 10 20 cm~ 2 (Stark et al. 1992). 

b H21 = (M — S)/{M + S) and H31 = (H — S)/(H + S) where S, M and H are the counts in a soft (0.3 - 1 keV), medium 
(1 — 2 keV), and hard (1 — 10 keV) energy bands respectively. 

1 Two optical candidates within matching radius. 

2 Elongated X-ray source located within three adjacent optical sources. 
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TABLE 2 

LMXB Luminosity Function Parameters 



Parameter 




Value 






All 


GCs 


Field 


Broken Power Law 








a\ 


— 1.10 ± U.zo 


1 io J. n n 


— l.oo ± U.zo 




-1.86 + 0.54 


-4.8+ 1.7 


-4.89 + 1.1 


l°g(L s [crgs i J) 


38.30 + 0.42 


38.49 ± 0.18 


38.34 + 0.15 


Truncated Power Law 








7 


-2.14 + 0.14 


-2.07 + 0.20 


-2.36 + 0.33 




TABLE 3 






Model Metallicities 


and Colors for Globular Clusters' 1 


[Fe/H] 


(S475 - Z850) 


(V - zg 50 ) 




(dex) 


(AB mag) 


(AB mag) 




-2.3 


0.851 


0.493 




-1.7 


0.907 


0.512 




-0.7 


1.243 


0.763 




-0.4 


1.425 


0.901 




+0.0 


1.588 


1.022 




+0.4 


1.897 


1.257 





Obtained from the Bruzual & Chariot (2003) models assuming an age of 13 Gyr 



